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Abstract 

The g-state Potts model partition function (equivalent to the Tutte polynomial) for a lat- 
tice strip of fixed width L y and arbitrary length L x has the form Z(G, q, v) = Yl!j=f' X c z,G,j{^z,G,jy 
where v is a temperature-dependent variable. The special case of the zero-temperature an- 
tiferromagnet (v = —1) is the chromatic polynomial P(G,q). Using coloring and trans- 
fer matrix methods, we give general formulas for Cx,g — Ylj=i G ' X c x,G,j f° r X = Z,P on 
cyclic and Mobius strip graphs of the square and triangular lattice. Combining these with 
a general expression for the (unique) coefficient cz,c,j of degree d in q: = U 2 d(^Y), 
where U n (x) is the Chebyshev polynomial of the second kind, we determine the num- 
ber of \z,g,/ s with coefficient in Z(G,q,v) for these cyclic strips of width L y to be 
nz(L y , d) = (2d + l)(L y + d+ for < d < L y and zero otherwise. For both cyclic 

and Mobius strips of these lattices, the total number of distinct eigenvalues Xz,gj is calcu- 
lated to be N Zi L y ,x = Cl v )- Results are also presented for the analogous numbers n P (L y ,d) 
and N P ^ Lyt \ for P(G,q). We find that n P (L y ,0) = n P (L y — 1,1) = Ml v -i (Motzkin num- 
ber), nz(Ly,0) = Cl v (the Catalan number), and give an exact expression for N PjLyt \. Our 
results for Nz,l v ,\ and Np t L y ,x apply for both the cyclic and Mobius strips of both the square 
and triangular lattices; we also point out the interesting relations Nz,l v ,x — 2N D A,tri,L y and 
Np,L y ,x = 2NDA,sq,Ly, where Nr>A,k,n denotes the number of directed lattice animals on the 

lattice A. We find the asymptotic growths N z .L y ,x ~ L y 1/2 ^ Ly and N PjLy> x ~ L y 1/2 3 Ly 
as L y — ► oo. Some general geometric identities for Potts model partition functions are also 
presented. 
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1 Introduction 



The g-state Potts model has served as a valuable model for the study of phase transitions 
and critical phenomena [1, 2]. On a lattice, or, more generally, on a (connected) graph G, 
at temperature T, this model is defined by the partition function 

Z(G,q,v) = J2^ H (1-1) 

with the (zero-field) Hamiltonian 

n = -jJ2^ (i.2) 

(ij) 

where a, = 1, q are the spin variables on each vertex i e G; (3 = {k B T)~ l ; and (ij) denotes 
pairs of adjacent vertices. The graph G = G(V, E) is defined by its vertex set V and its edge 
set E; we denote the number of vertices of G as n = n(G) — \V\ and the number of edges of 
G as e(G) = \E\. We use the notation 

K = PJ , a = u" 1 = e K , v = a-l (1.3) 

so that the physical ranges are (i) a > 1, i.e., v > corresponding to oo > T > for the 
Potts ferromagnet, and (ii) < a < 1, i.e., — 1 < v < 0, corresponding to < T < oo for 
the Potts antiferro magnet. One defines the (reduced) free energy per site / = —f3F, where 
F is the actual free energy, via 

f({G},q,v)= hm \n[Z(G,q,vy/ n ] . (1.4) 

where we use the symbol {G} to denote lim^oo G for a given family of graphs. 

Let G' = (V, E') be a spanning subgraph of G, i.e. a subgraph having the same vertex 
set V and an edge set E' C E 1 . Then g, v) can be written as the sum [3]- [6] 

Z(G,q,v) = Y.^'^ 
G'CG 



n{G) e(G) 
r=fc(G) s=0 

where k(G') denotes the number of connected components of G' and z rs > 0. Since we only 
consider connected graphs G, we have k(G) = 1. The formula (1.5) shows that Z(G,q,v) 
is a polynomial in q and t> and enables one to generalize q from Z + to M + and, indeed, to 
C. The Potts model partition function Z(G,q,v) on a graph G is essentially equivalent to 
the Tutte polynomial [7]- [9] T(G,x,y) and Whitney rank polynomial R(G,£,r)) [2], [4], [10] 
for this graph, as discussed in the appendix. 
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One special case of the Potts model partition function that is of particular interest is the 
zero-temperature limit of the Potts antiferromagnet (AF). For sufficiently large q, on a given 
lattice or graph G, this exhibits nonzero ground state entropy (without frustration). This is 
equivalent to a ground state degeneracy per site (vertex), W > 1, since S = k B \nW. The 
T = (i.e., v = —1) partition function of the above-mentioned 5- state Potts antiferromagnet 
(PAF) on G satisfies 

Z(G,q,-l)=P(G,q) (1.6) 

where P(G, q) is the chromatic polynomial (in q) expressing the number of ways of coloring 
the vertices of the graph G with q colors such that no two adjacent vertices have the same 
color [3, 10, 12, 13]. This is termed a proper vertex coloring of G. The minimum number of 
colors necessary for this coloring is the chromatic number of G, denoted x(G). Thus 1 

W({G},q) = lim P{G,q) 1/n . (1.7) 

A second special case of the Potts model partition function is for infinite temperature, i.e., 
K = or equivalently v — 0. In this case, as is clear from either (1.1) with (1.2) or from 
(1.5), Z reduces to the single term 

Z(G,q,0) = q n ^ . (1.8) 

A third special case is the Potts ferromagnet in the limit of zero temperature, T — ■> 0, i.e., 
v — > 00. 

Here we shall consider families of graphs that are cyclic strips of the square and triangular 
lattices, taken to be oriented horizontally, with free transverse boundary conditions, denoted 
FBCy, and periodic longitudinal boundary conditions, denoted PBC X . A given strip of the 
triangular lattice may be visualized by starting with the corresponding strip of the square 
lattice and adding diagonal bonds joining, say, the upper left to lower right vertices of each 
square. The length and width are taken to be L x and L y vertices, and the total number of 
vertices is n = L x L y . We shall also make some comments about the analogous strips with 
twisted periodic longitudinal, i.e., Mobius, boundary conditions, denoted TPBC X or Mb. 
These various families of lattice strip graphs are examples of recursive families of graphs, in 
the sense that a strip of length L x + 1 is constructed by the addition of a given subgraph 
(here, a transverse layer of the strip) to the strip of length L x . As derived in [15, 16], using a 
type of transfer matrix argument, a general form for the Potts model partition function for 
the strip graphs considered here, or more generally, for a recursively defined graph comprised 
of L x = m repeated sub units (transverse layers of the strip here), is 

Z(G,q,v)= Cz,GAXz,G,) m (1.9) 

1 At certain special points q s (typically q s — 0, 1, .., x(G)), one has the noncommutativity of limits 
lim q ^ qs linin^oo P(G, q) 1 /™ =/= \im n ^ 00 lim q ^ qs P(G, q) 1 ^ 1 , and hence it is necessary to specify the order 
of the limits in the definition of W({G}, q s ) [14]. 
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where Xz,G,j is a function of q and v, and both Xz,gj an d the coefficients cz,gj are independent 
of L x . The coefficients cz t G,j are functions only of q and, indeed, are polynomials, for the 
cyclic strip graphs of interest here. (We shall comment on other types of strip graphs below.) 
From (1.6), it also follows that the chromatic polynomial has the same structure, 

P(G,q)= CP,oA*P,Oj) m (1-10) 

J'=l 

where Np t c,\ depends on G; cp t a,j and Xp,G,j are independent of L x = m; and 

N p , g ,x<N z ,g,x . (1.11) 

As will be shown below, for the strips under consideration here, this inequality is realized 
as an equality for the circuit graph C n and as a strict inequality for widths L y > 2 (see 
Tables 1 and 3). For the same reason, it follows that when one sets v — — 1 in Z(G,q,v), 
i.e. specializes to the T = Potts antiferromagnet, a certain number Nz,g,x ~ ^p,g,a °f the 
Azg/s vanish and the remaining, nonvanishing XzGj(Q, V T S are precisely those occurring in 
P(G,q): 

If \ z ,Gj(q,v = -1) ^ then A Z)G j(<Z, -1) = A PlGJ (g) . (1.12) 

Previous calculations of chromatic polynomials for recursive families of graphs of arbitrary 
length include [18]- [49]; we shall also use the results of calculations of Potts model partition 
functions (Tutte polynomials) for cyclic and Mobius strip graphs of arbitrary length in 
[15, 16, 47]. Some relevant results on transfer matrices and the Temperley-Lieb algebra are 
in [29], [50]-[53]. 

Several basic questions about the structure of Z(G, q, v) and P(G, q) for these cyclic strip 
graphs are the following: 

1. Can one obtain a general formula for the coefficients Cz,G,j an d Cp,G,j m (1-9) and 
(1.10)? 

2. What are the respective sums of the coefficients c Z) c,j a nd cp t c,j in (1.9) and (1.10)? 

3. For a given cyclic strip G, how many Xz,g/ s m (1-9) have a particular coefficient, and 
how many X P: g/s in (1.10) have this coefficient? 

4. What are the total numbers of different terms N z ,g,x a nd N PiG ,x in (1.9) and (1.10)? 

In this paper, we shall obtain a general formula for the coefficients cz,G,j an d cp : G,j, a n 
answer to the second pair of questions, and, based on these results, answers to the remaining 
two pairs of questions. We shall also comment on other types of lattice strips. 

The chromatic polynomial, can be calculated by various methods, including iterative 
application of the deletion-contraction theorem (e.g., [10]), a certain matrix technique [21], 



3 



a generating function method [33, 40], or a coloring (compatibility) matrix method [22, 
43]. Similarly, the Potts model partition function or equivalent Tutte polynomial can be 
calculated by iterative application of the generalized deletion-contraction theorem [10] or by 
transfer matrix methods including a relation with the Temperley-Lieb algebra [29, 51, 53]. 
Of course these different methods can be used to provide cross-checks; for example, iterated 
deletion-contraction and transfer matrix methods were used to check each other in [16, 47] 
and deletion-contraction and coloring matrix methods were used together in works such as 
[45]. The last, the coloring matrix method, is also useful for deriving rigorous upper and 
lower bounds on q) [22], [23], [24]. In this method, one first selects a transverse slice 

of the strip, denoted Cl v - For our cyclic graphs, this is simply a line graph with L y vertices 
(oriented vertically, given that we orient the long direction of the strip horizontally). Denote 
an allowed g-coloring of this path as c(Cl v )- The number of allowed colorings of the path 
Cl v is M = P(£L y ,q)- For our chromatic polynomials of cyclic strips of the square and 
triangular lattices, P(C n , q) = P(T Ly , q), where T n is the tree graph with n vertices and 

P(T n ,q) = q(q-l) n - 1 . (1.13) 

Now focus on two adjacent paths Cl v and C' Ly . Define compatible g-colorings of these 
paths as colorings such that no two adjacent vertices v G Cl v and v' G C' Ly (i.e. vertices 
connected by an edge = bond of the lattice strip graph) have the same color. One can then 
associate with this pair of paths an M x M dimensional symmetric matrix T with entries 
%(c Ly ),c(C' L ) = 1 or if the g-colorings of Cl v and C' Ly are or are not compatible, respectively. 
Then the chromatic polynomial of the cyclic strip of the lattice A, taken here to be square 
(sq) or triangular (tri), is (with L x = m) given by 

P(A, L y x L x , FBCy, PBC X , q) = Tr(T m ) . (1.14) 

Since T is a symmetric real matrix (indeed, composed only of 0's and l's, although we do not 
need this here), it can be diagonalized by an orthogonal transformation, so that the above 
trace is Tr(T m ) = Yl,j(^P,G,j) m -, where the Ap iG j's are the eigenvalues of T. In the context 
of the g-coloring problem, if one denotes the multiplicity of the j'th distinct eigenvalue 
by Cp t c,j, one obtains the formula (1.10). That is, from the coloring matrix viewpoint, 
c p,gj is the dimension of the invariant subspace in the full A/"-dimensional space of coloring 
configurations of the transverse slice of the strip corresponding to the eigenvalue \p,c,j- For 
the full temperature-dependent Potts model partition function, one can define an analogous 
matrix, T z whose entries, rather than being and 1, are appropriate Boltzmann weights 
for the given spin configurations on the successive transverse slices Cl v and C! Ly . Then this 
partition function is (with L x = m) 

Z(A,L y x L x ,FBC y ,PBC x ,q,v) = Tr[(T z ) m ] (1.15) 

which was used in [15, 16] to derive the formula (1.9). Note that in this case all colorings, 
including those that yield the same color on adjacent vertices, are allowed. One can obtain 
these formulas (1.14) and (1.15) for sufficiently large integral q that the multiplicities cp jGj - 
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and cz,gj are positive integers; in both of these cases, the fact that the coefficients are 
multiplicities shows that for sufficiently large positive integer q they are also positive integers. 
In the case of the chromatic polynomial, it is obvious that cp t a,j can depend only on q; in 
the case of the full Potts model partition function, for the cyclic strip graphs considered 
here, this also follows from (1.15) and (1.9), since the multiplicities of the eigenvalues cannot 
depend on the variable parameter v G (—1, oo). 

Note that, while the coefficients cp t a,j and cz,G,j can be obtained as multiplicities of 
the distinct eigenvalues \p,c,j and Xz,gj f° r sufficiently large integer q, when one considers 
positive q < 4, they may be zero or negative (see eqs. (2.17)-(2.21) below). More generally, 
while they play the role of eigenvalue multiplicities for sufficiently large integer q, the domain 
of their definition may be generalized via (1.5) to q G M+ or, indeed, to q G C. Some 
properties of determinants of coloring matrices for various strip graphs are given in the 
appendix. 

The dimension of the space of coloring configurations, A/", is equal to the sum of the 
multiplicities of each distinct eigenvalue, i.e., the sum of the dimensions of the invariant 
subspaces corresponding to each of these distinct eigenvalues. For the chromatic polynomial, 
this is Af, which is equal to the sum 

Cp, G = C P0j (1-16) 

i=i 

while for the full Potts model partition function, we shall denote it as 

C Z , G = c *Ai • (1-17) 

3=1 

General results for these sums will be given below for the strips of interest. 

Before proceeding, we note that for the chromatic polynomial, since two identical color 
configurations on Cl v and C ' l v are incompatible, the diagonal elements of the coloring matrix 
are zero and hence its trace is zero: 

Tr(T(G,q)) = for P(G,q) . (1.18) 

Since the coefficients cp t Q,j are just the multiplicities of the given eigenvalues Xp,G,j, we can 
write this as the sum over distinct Ap ; gj's: 

'■'■■<;,<\r.<;, = . (1.19) 

3=1 

In the following, since we shall present results that are applicable to cyclic strips of both 
the square and triangular lattices, we shall often leave the lattice type A and the boundary 
conditions (FBC y , PBC X ) implicit in the notation where they are obvious. The fact that 
these results will be applicable to cyclic strips of both the square and triangular lattices 



5 



depends on the property that (having expressed the strip of the triangular lattice as a 
strip of the square lattice with diagonal bonds added to each square, as described above) 
the transverse slices, i.e., line graphs of L y vertices, of both the cyclic square-lattice and 
triangular-lattice strips are identical. 2 Similarly, for quantities that are independent of some 
part of G, such as cz,G,j and ^z,gj, which are independent of L x , this will be incorporated 
in the notation. 



2 Coefficients in Potts Model Partition Functions for 
Cyclic Lattice Strips 

In this section we address the first pair of questions posed in the Introduction. From our 
exact solutions of the chromatic polynomials and, more generally, of the full temperature- 
dependent Potts model partition functions on cyclic strips of the square and triangular 
lattice, we have found that the coefficients cz,G,j and cp t c,j are polynomials in q. Further, we 
have found that for a given cyclic strip of the square or triangular lattice with width L y , the 
coefficients cp t c,j and cz,G,j are of a limited set; only one type of polynomial of each degree in 
q occurs 3 , and, denoting this as c^ d \ there are coefficients of the form with < d < L y . 
We infer the following general formula for c^ d \ i.e. for the multiplicity of the corresponding 
eigenvalue of T z : 

where U n (x) is the Chebyshev polynomial of the second kind, defined by (e.g. [54, 55]) 

U n {x) = Ira 



[x + iVT^x 1 ) n+1 ] 



3=0 



_ {2x) n ~ 13 (2.2) 



where in eq. (2.2) and similar equations below, the notation [|] in the upper limit on the 
summand means the integral part of |. In the first line of (2.2), x must be real, but the 

2 Although the transverse slices are identical for the cyclic strips of the square and triangular lattice, they 
are different for other lattices, such as the honeycomb = brick lattice, and this difference was also significant 
for the earlier use of coloring matrix methods to obtain rigorous bounds on the Potts model ground state 
degeneracy per site, W [2 2] -[24]. 

3 From the studies of Mobius strips of the square lattice [18, 42, 44, 48], we find that the coefficients are 
again polynomials, of the form ±c( d '; however, for Mobius strips of the triangular lattice, the coefficients 
are not, in general, polynomials, but algebraic functions of q, as was shown by the exact solutions for the 
chromatic polynomial in [42] and for the full Potts model partition in [47] for the L y — 2 Mobius strip. For 
strips with torus boundary conditions, exact solutions [45, 46] show that the coefficients are polynomials of 
degree d in q, but it is not, in general, the case that there is a unique coefficient of degree d and they are 
not, in general, of the form c^ d \ We shall comment further on this below. 
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polynomial in the second line of (2.2) constitutes a definition of U n (x) for complex x. The 
first few of these coefficients are, with x — 

c (0) = U (x) = 1 (2.3) 

c « = C/ 2 ( x ) = 4a; 2 - 1 

= 5-1 (2-4) 

c (2) = C/ 4 (a;) = 16a; 4 - 12a; 2 + 1 

= g 2 -3g+l (2.5) 

c ( 3 ) = u 6 (x) = 64a; 6 - 80a; 4 + 24a; 2 - 1 

= g 3 - 5q 2 + 6q - 1 (2.6) 

c ( 4 ) = f/ 8 ( x ) = 256a; 8 - 448a; 6 + 240a; 4 - 40a; 2 + 1 
= q A - 7q 3 + I5q 2 - lOq + 1 

= (g-l)(g 3 -6g 2 + 9g-l) (2.7) 
and so forth for higher d. We have found, as an equivalent formula [46, 48] 

c {d) = f[(q - q d , k ) (2.8) 
k=i 

where 



,^2 + 2cos(^)=W(^), for *=l,2,...d. (2.9) 

Note that the apparent square root in (2.1) is absent in the actual formula (2.1) for the 
c^'s, which applies for arbitrary complex q. One way to derive (2.1) for cyclic strips is to 
use the fact that the partition function is given by the trace (1.15), which, with (1.9) makes 
clear the role of the coefficients cz,gj as multiplicities of eigenvalues of the transfer matrix 
Tz; then, one utilizes the connection with the Temperley-Lieb algebra to infer that these 
multiplicities are given by (for sufficiently large integer q where the are positive 
integers), and finally, one continues this result to general complex q. Related discussions 
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of dimensions of invariant subspaces of operators in the Temperley-Lieb algebra are given 
in [52], [53]. The form (2.8) is also related to the property that each coefficient has a 
simple zero at the Tutte-Beraha number B 2d+ i] indeed, one could derive the form (2.1) by 
using this as a starting point. In this approach, one would start with the property that 
has the factor (q — B 2d+ i) — (q — Qd,i)- Since B 2d+ \ is, in general, irrational, it is 
necessary to symmetrize the product = Ylil ~ Qz) ° ver the zeros q z in order to obtain 
rational (indeed, integer) coefficients for the c^ d \ This symmetrization yields (2.8). To 
see this, we note that to incorporate this zero of at q — B 2 d+i one can start with the 
identity [exp(2?™/(2d + 1))] M+1 -1 = 0, express this as [cos(^) + % sin(^r[)} 2d+1 , expand 
the resulting expression, and then impose the condition that the polynomial vanishes at 
Q — Qd,i — B 2 d+i by setting cos(27r/(2<i + 1)) = q/2 — 1. This actually yields two equations, 
one for the real part and one for the imaginary part. Concentrating on the real part, we 
observe that the resultant equation has the form (l/2)(g — 4)(c^) 2 . This yields the result 



c (d) 



— (TWi(f-l)-l)] (2.10) 
where T n (x) is the Chebyshev polynomial of the first kind, defined by (e.g., [54, 55]) 

T n (x) = Re[(x + iVl -x 2 ) n 

[-] 

As with U n (x), in the first line of (2.11), x must be real, but the polynomial in the second 
line of (2.11) constitutes a definition of T n (x) for complex x; note also that the right-hand 
side of the second line is defined to be 1 for n — 0. Simplifying eq. (2.10), one finally obtains 
eq. (2.1) in terms of the Chebyshev polynomial of the second kind. The starting point is 
equivalent to the identity [exp(27rk/(2d + l))] 2d+1 — 1 = 0, and hence it follows that 
vanishes at qa : k for 2 < k < d as well as for k — 1. This implies the formula (2.8). 

The general formula (2.1) or the equivalent result (2.8) gives a deeper understanding 
of, the known exact solutions for chromatic polynomials and the full Potts model partition 
functions on cyclic strips of the square and triangular lattice. For reference, these include 
the degenerate case L y — 1 (circuit graph), and the cases L y = 2 in [18], L y = 3 in [42], 
and L y = 4 in [48] for cyclic strips of the square lattice; and the cases L y = 2 in [42] and 
L y = 3 and L y = 4 in [46] for cyclic strips of the triangular lattice. For the full Potts model 
partition functions, the exact solutions include the elementary L y — 1 solution for the circuit 
graph, and the L y = 2 solutions in [16] and [47] for cyclic strips of the square and triangular 
lattices, respectively. Explicit examples will be given below. As noted, the coefficients with 
degree d occur for < d < L y . In the following we shall work out implications of the formula 
(2.1) for the structure of the chromatic polynomial and full Potts model partition function 
for arbitrarily wide cyclic strips of the square and triangular lattices. 
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From the formula (2.1) or (2.8) and properties of the Chebyshev polynomials of the 
second kind, one can establish a number of useful properties of the coefficients c^ d \ First, 
we recall that a generating function for the U n (x) is (e.g. [54]) 

71=0 

The formal sum in (2.12) converges if \x\ < 1 and \z\ < 1; however, for our purposes, we 
shall only use it as a means of extracting the coefficient of the z n term as U n (x) ) and this 
can be done for any x and z. We infer the following properties of the c^: 

1. is a polynomial of degree d in the variable q, with integer coefficients, highest-order 
term equal to q d , and alternating signs for the subsequent terms of descending degree 
in q. 

2. We recall the definition of a unimodal polynomial as one with the property that the 
magnitudes of its coefficients increase monotonically up to a point and then (with 
possible equality of the maximal two coefficients) decrease monotonically. From the 
property that the magnitudes of the coefficients of U n (x) are unimodal, it follows 
directly that the same is true of the magnitudes of the coefficients of c^ d \ 

3. The zeros of are real, positive, simple, lie in the interval < q < 4 and are given 
as follows: 

C W=0 at q = q d , k , k = l,2,...,d (2.13) 
In the limit d — > oo, these zeros become dense in this interval. 

4. If d > 1, then one of the zeros of is at a Tutte-Beraha number, namely, 

q d ,i = B 2d+1 (2.14) 
where the Tutte-Beraha number B r is defined as [9], [26] 

5 r = 2 + 2cos(y) =4cos 2 ^) ,r = 1,2,... (2.15) 

The first Tutte-Beraha numbers are B\ = 4, B 2 = 0, -B3 — 1, B4 — 2, B 5 = (3 + 
VE )/2 = 2.618.., B % = 3. 

5. Depending on the value of d, may have other zero(s) at Tutte-Beraha number(s) 
distinct from B 2d +i in (2.14). This cannot happen if 2d + 1 is prime. If 2d + 1 is not 
prime, consider the case where there exists d' such that k(2d' + 1) = 2d + 1. Since 
2d + 1 and 2d' + 1 are both odd, it follows that k is also odd, so we can write this 
equation in a symmetric manner as 

(2d' + l)(2rf" + l) = 2d+l . (2.16) 
9 



(2.12) 



Then has additional zero(s) at B 2 d>+i and B 2 d"+i- These degenerate into a single 
additional zero if 2d + 1 = p 2 where p is prime. Note that the factorization (2.16) 
is not unique. For example, consider for d = 52. Then 2d + 1 has the full and 
unique factorization 105 = 3-5-7. There are thus three ways of writing the twofold 
factorization in (2.16): 3 • 35, 5 • 21, and 7 ■ 15, and hence, in addition to the zero at 
the Tutte-Beraha number B 2 d+i = B W5 , the coefficient c 1 - 105 ^ also has zeros at the six 
other Tutte-Beraha numbers B 3 , B 5 , B 7 , Z? 15 , B 2 \, and -B35. As this example makes 
clear, in other cases there could be more than three twofold factorizations of the form 
(2.16). 

6. As a particular case of the previous item, if d — 1 mod 3, then, setting d — 3j + 1 and 
substituting in the above equation, we find (2d' + I) (2d" + 1) = 3(2 j + 1). If j — so 
that d = 1, i.e., = q — 1, then there is only the single zero at B 3 = 1 given by (2.14). 
If j = 1, i.e., d = 4, then has a zero at B 3 and at B 9 . If j > 2, then in addition to 
the zero at i?2d+i given by (2.14), c {d) has zeros at the Tutte-Beraha numbers -B3 and 
B 2 j+\ (as well as possible others). Note that since has only simple zeros, it follows 
that if d — 1 mod 3, then has a simple factor (q — 1). (The expression given above for 

illustrates this.) 

7. Since U 2n (0) = (-l) n , it follows that 

c (d) = for g = Q ( 217 ) 



If q = l then c (d) = <j if d = 1 mod 3 (2.18) 



9. 



10. 



for 


q = 






ifd = 


= mod 3 


{"., 


ifd = 


= 1 mod 3 




if d = 


= 2 mod 3 




if d = 


0, 1 mod 4 




if d = 


2,3 mod 4 


1 


if d = 


0,2 mod 6 




if d = 


1 mod 6 




if d = 


3, 5 mod 6 




if d = 


4 mod 6 



If , = 2 then C «= ' "I ' PIS) 



11. 



If q = 3 then = 2 , • , = o K , a (2-20) 



c ( rf ) = 2rf + 1 for q = 4 . (2.21) 



12. As a consequence of the inequality |C4(x)| < n + 1 for — 1 < x < 1, we have the 
inequality 

\c (d) \<2d+l for0<g<4. (2.22) 
This inequality is saturated at q — 4, i.e. x = 1. 
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13. Using the recursion relation for Chebyshev polynomials, 

U n+1 (x) = 2xU n (x) - U n -!(x) (2.23) 

iteratively, we find the resulting recursion for the c^ d \ 

c (d+i) = ( q _ 2 y{d) _ c (d-i) _ ^ 2 2 ^ 



If one lets 

/ 9\ 

g = 2 + 2cos# = 4cos 2 (^-J (2.25) 

one sees that the argument of the Chebyshev polynomial of the second kind in (2.1) is given 
hy x = ^ = cos(#/2). A relevant identity is (with uj = 9/2 here) 

. . sin((n + l)u>) 

UJcosuj) = — ^ '—L . 2.26 

sincj 

The fact that the transformation (2.25) applies, with 9 real, for q G [0,4] shows the special 
role of this interval for Potts model partition functions on the strips under consideration 
here. This makes an intriguing connection with the 2D Potts ferromagnet, which has a 
second-order phase transition for the values q = 2 (Ising), q — 3, and q — 4; and a first-order 
phase transition for q > 5. Of course, the physical thermodynamic behavior of the Potts 
ferromagnet on the infinite-length, finite width strips is quite different from that of the model 
on 2D lattices; in the former cases, since these are quasi-one-dimensional, the model has only 
a zero-temperature critical point rather than a finite-temperature phase transition. For the 
antiferromagnetic Potts model, the transformation means that the value q = 4 is special (as 
the positive value of q where the angle 9 changes from being real to imaginary) and is in 
accord with the special role of q = 4 in this model on the infinite two-dimensional lattice 
[29] (see also [32]). This is somewhat similar to earlier situations in which studies of lower- 
dimensional realizations of spin models gave insight into properties of higher-dimensional 
realizations (e.g. [56]- [58] for O(N) models and [16, 47, 59] for Potts/Ising models). However, 
in assessing the connection of the special role of the value q = 4 in the coefficients with 
the T = critical point at q = 4 for the Potts antiferromagnet on the triangular lattice, one 
must also taken into account the fact that the formula (2.1) also applies for the cyclic strips 
of the square lattice that have been studied; however, again taking the limit L y — > oo, the 
Potts antiferromagnet on the square lattice has a zero-temperature critical point at q = 3 
rather than at q = 4 [17]. 

One of the interesting aspects of the zeros q^k of is that at these values of q, the 
critical exponents of the 2D Potts model are rational. We recall that the thermal, magnetic, 
and tricritical exponents for the paramagnetic to ferromagnetic transition in the 2D g-state 
Potts ferromagnet are, for q < 4, [60] -[62], [2] 

y* = (2 ' 27) 



n 



(3 — it) (5 — u) , n 

y h , y h ,tnc. = 1 A{ 2 _u) ( } 

where 

u = -arccosf^) (2.29) 

7T V 2 / 

so that < it < 1 for 4 > q > 0. The angle 6 in (2.25) is equal to nu in (2.29), and hence 

2k 

q = q d ,k <£> u = ^—^ . (2.30) 

Note that the converse does not hold; that is, there are rational values of u that are not of 
the form (2.30) and hence do not correspond to any of the q^k- An example is u — 1/2. 



3 Determination of np(L y ,d) for Cyclic Strips of the 
Square and Triangular Lattices 

In this section and the next we use (2.1) together with two theorems (eqs. (3.10 and (4.1) 
below) to determine structural properties of the chromatic polynomial and the full Potts 
model partition function for cyclic strip graphs G of the square and triangular lattices. Let 
us define np(L y ,d) as the number of terms Xp,G,j in P(G,q) that have as their coefficients 
°p,G,j — an d nz(L y ,d) as the number of terms Xz,G,j i n Z(G,q,v) that have as their 
coefficients cz,gj — c<yd ^ ■ For the cyclic strip graphs under consideration here, these coeffi- 
cients are independent of L x and depend on L y and d; furthermore, they are the same for 
both square and triangular strips, as discussed further below. While the individual Ap^j's 
in (1.10) are, in general, different for the cyclic strips of the square and triangular lattices, 
the total number of \p,g/ s is the same, so we use the short notation N PjLy ^. The same is 
true for the full partition function, so we use the notation Nz,L y ,\- 

Since the chromatic polynomial P(G, q) is the special case v — — 1 (i.e. zero-temperature 
antiferromagnet) of the general (finite-temperature, J positive or negative) Potts model 
partition function, (1.6), it follows that 

n P (L y ,d) < n z (L y ,d) . (3.1) 

The total number, Np t c,x, of different terms Xp,G,j in (1-10) is given by 

Ly 

N P , Ly , x = J2ML y ,d) . (3.2) 

d=0 

Since each term Xp,G,j is a distinct eigenvalue of the coloring matrix T, the number iVp^ ^ is 
the number of different invariant subspaces in the full A/-dimensional space of coloring con- 
figurations of the transverse slices of the strips. Similarly, for the full Potts model partition 
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function, 

Ly 

N z , Ly , x = J2ML y ,d) . (3.3) 

d=0 

and analogously this represents the number of different invariant subspaces for the matrix 
Tz- It was shown in [44] that the \p t G(L y ),j's are the same for the cyclic and Mobius strips 
(although the corresponding cp^j's are different) 

^P,G(L y ),FBCy,PBC x ,j = ^P,G(L y ),FBCy,TPBC x ,j V j (3.4) 

and hence the total number of terms is also the same: 

Np t G(L y ),FBCy,PBC x ,\ = Np t G(L y ),FBCy,TPBC x ,\ ■ (3.5) 

The argument in [44] relied upon the local nature of the deletion-contraction operations used 
in calculating P(G, q), and the same property is true of the deletion-contraction theorem used 
for calculating the Tutte polynomial, or equivalently, the Potts model partition function, so 
that the following generalizations of (3.5) holds [15, 16]: 

^Z,G(L y ),FBCy,PBC x ,j = ^Z,G(L y ),FBCy,TPBC x ,j V j (3.6) 

and hence 

Nz,G{L y ),FBC y ,PBC x ,\ = N Z ,G(L y ),FBCy,TPBC x ,\ ■ (3.7) 

For the sum of the coefficients in (1.10), i.e., (1.16), we have 

C P , Ly = cp^ = 5> P (L tf ,d)cM (3.8) 

j=l d=0 

and for the corresponding sum of coefficients in (1.9), 

Nz,Ly,\ Ly 

C z , Ly = cz,L v j = Y,n z {L v ,d)cV> . (3.9) 

3=1 d=0 

We first recall a theorem specifying Cp^ y for cyclic strips of the square and triangular 
lattice [44]: 



Theorem 1. 

Cp,L y = P(T Ly ,q)=q(q-l) L y- 1 . (3.10) 

Proof. Using coloring matrix methods [22, 43], one has that Cp t L is equal to the 
chromatic polynomial for the coloring of the transverse slice of the strip, which is a line 
with L y vertices. This is a special case of a tree graph, for which the elementary general 
result in eq. (1.13) holds. □ 
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Next, we have 



Theorem 2. The np(L y , d), d = 0, 1, ..L y are determined as follows. One has 

n P {L y ,d) = for d > L y , (3.11) 
np(L y ,L y ) = l (3.12) 

and 

ra P (l,0) = l (3.13) 
with all other numbers np(L y ,d) being determined by the two recursion relations 

n P (L y + 1,0) = n P (L y ,l) (3.14) 

and 

np(L y + l,d) = np(Ly,d—l)+np(Ly,d)+np(Ly,d+l) for L y > 1 and 1 < d < L y +1 . 

(3.15) 

Proof. We substitute for from eq. (2.1) in eq. (3.10). We obtain another equation by 
differentiating this with respect to q once; another by differentiating twice, and so forth up 
to L y -fold differentiations. This yields L y + 1 linear equations in the L y + 1 unknowns, 
n P (Ly, d), d — 0, 1, L y . We solve this set of equations to get the n P (L y , d). □ 

A corollary is that 

n P (Ly,L y -1) = L y . (3. 16) 

The numbers np(L y ,d) can be viewed as integer sequences in L y for a given value of 
d. We find that (for L y > 1 where our strips are defined) np(L y ,0) is a Motzkin number 



np(L y ,0) = M Ly _ 1 (3.17) 
where the Motzkin number M n is given by 

M n = J2(-iyC n+1 -J n ) (3.18) 

3=0 ^ ^ 

where 

a, = 4^") (3.i9) 

n + 1 \ n J 

is the Catalan number. (We also use the symbol C n to denote the circuit graph with n 
vertices; the meaning will be clear from context.) The Catalan and Motzkin numbers oc- 
cur in many combinatoric applications [63]- [69]. Among these is the construction of non- 
intersecting chords on a circle; the number of ways of connecting a subset of n points on a 
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Table 1: Table of numbers np(L y ,d) and their sums, Np t L v ,x f° r cyclic strips of the square and triangular 
lattices. Blank entries are zero. 



T 1 7 

L y I d — ► 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Np, Ly ,X 


1 


1 


1 




















2 


2 


1 


2 


1 


















4 


3 


2 


4 


3 


1 
















10 


4 


4 


9 


8 


4 


1 














26 


5 


9 


21 


21 


13 


5 


1 












70 


6 


21 


51 


55 


39 


19 


6 


1 










192 


7 


51 


127 


145 


113 


64 


26 


7 


1 








534 


8 


127 


323 


385 


322 


203 


97 


34 


8 


1 






1500 


9 


323 


835 


1030 


910 


622 


334 


139 


43 


9 


1 




4246 


10 


835 


2188 


2775 


2562 


1866 


1095 


516 


191 


53 


10 


1 


12092 



circle by non-intersecting chords is M n , while the number of ways of completely connecting 
2n points on the circle by such chords is C n . Summing over subsets of points connected by 
chords, this yields the well-known relation 




(3.20) 



As a consequence of our relation (3.14), the equation (3.17) also implies that 

n P (L y , 1) = M Ly . (3.21) 

A generating function is 

G np{Lyfi) {x) = i-s-Vij 2 *-^ = £ np{L y, O)^- 1 . (3.22) 

Ly=l 

In eq. (5.1.3) below we shall give an exact determination of the total number, Np t L \, of 
types of terms Xp t L y ,j- In Table 1 we list the numbers np(L y , d) and the sums N PjLyt \ for the 
first several widths of cyclic strips of the square and triangular lattice, 1 < L y < 10. Others 
can easily be calculated using our general formulas. 

Now let us consider a random walk on the nonnegative integers such that in each step 
the walker moves by +1, —1, or units. Denoting m(n, k) as the number of walks of length 
n steps starting at and ending at k, we obtain the Motzkin triangle in Table 2. The first 
column, corresponding to k = is the number of walks defined above that return to the 
origin after n steps. This is given by the Motzkin number defined in (3.18); m(n, 0) = M n . 
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Table 2: Table of numbers m(n, k) for the number of a random walk from to k in n steps and the row 
sums S n . Blank entries are zero. 



n I k — > 





1 


2 


3 


4 


5 


6 


7 


8 


9 


<Sn 





1 




















1 


1 


1 


1 


















2 


2 


2 


2 


1 
















5 


3 


4 


5 


3 


1 














13 


4 


9 


12 


9 


4 


1 












35 


5 


21 


30 


25 


14 


5 


1 










96 


6 


51 


76 


69 


44 


20 


6 


1 








267 


7 


127 


196 


189 


133 


70 


27 


7 


1 






750 


8 


323 


512 


518 


392 


230 


104 


35 


8 


1 




2123 


9 


835 


1353 


1422 


1140 


726 


369 


147 


44 


9 


1 


4246 



The second column, m(n, 1), is given by the first differences of Motzkin numbers: 
The row sums in Table 2 will be important below; we denote them as 

n 

S n = ^m(n,k) . (3.24) 

k=0 

Notice that Table 1 can be viewed as the combination of two Motzkin triangles, i.e. Table 
2, as follows, 

rip(L y , d) = m{L y — 1, d) + m(L y — 1, d — 1) . (3.25) 

In light of the property that the coefficients (= multiplicities of eigenvalues of the transfer 
matrix) are Chebyshev polynomials of the second kind for the cyclic strips of the square 
and triangular lattices and for Mobius strips of the square lattice (discussed below), and the 
finding in eq. (3.17) for the numbers of Ap^j's with coefficients and it is interesting 
to note that in [65] a connection has been shown between what is called a Motzkin polynomial 
and this Chebyshev polynomial. To see this, one defines s^o — M h , s ft)1 = M h — M h _ x and 
shows that, with Sh >n = Sft,n-i _ Sft-i,n-i — s h~2,n-2 for 1 < u < h, it is possible to express 
Sh, n as Sh, n — a nMh + a n _\M h _\ + ... + a M h _ n , where the coefficients aj are independent of 
h. The Motzkin polynomial is defined as S n (x) = a n x n + ... + a±x + a , and the connection 
is that S n (x) = U n {^) [65]. 

Certain (g-independent) relations between the np(L y , d) can be derived by evaluating the 
sum (3.8) and its result (3.10) for q — 0, 1, and 2. Setting q = in (3.8) and (3.10), and 
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using (2.17), we have 



J2n P (Ly,d)c^(q = 0) = £(-1) d n P (L y , d) = . (3.26) 



d=0 d=0 

Next we evaluate (3.8) and (3.10) for q — 1. Yi L y = 1, this only involves one term, 
np(l, 0) = 1. So consider L y > 2 and set L y = 3k + r with integer A; > and r = 0, 1, or 
2. To express the equation compactly, it is convenient to use a Heaviside step function 9(z), 
where 9(z) — 1 if z > and #(2:) = if z < 0. Using (2.18), we have 

L y k k-l+9(r-l) 

J2n P (L y ,d)cW(q = l) = J2ML y ,3j)- J] n P (L y ,3j + 2) 

d=0 j=0 j=0 

= (3.27) 

where here and in other summations, it is understood that if the upper limit on the sum- 
mation is negative, the sum is zero. A similar relation can be obtained in a straightforward 
manner by substituting q = 2 in (3.8) and (3.10). Other similar relations among the np(L y , d) 
can be obtained by evaluating (3.8) and (3.10) for q = 3 and q = 4; for these relations, the 
right-hand sides grow with L y , in contrast to (3.26)-(3.27), where the right-hand sides are 
constant. 



4 Determination of nz(L yi d) for Cyclic Strips of the 
Square and Triangular Lattices 

We next generalize these results to the full Potts model partition function. We determine 
the numbers riz(L y ,d) by observing that in the cyclic strips under consideration here, the 
only difference between the coloring for the T = Potts antiferromagnet and for the general 
Potts model at finite temperature is the relaxation of the constraint that no two adjacent 
vertices can have the same color (indeed, for ferromagnetic coupling, there is a preference 
for these vertices to have the same color). Hence, one simply counts all possible colorings, 
and obtains 

Theorem 3. 

Ly 

J2nz(L y ,d)c^ = q L * . (4.1) 

d=0 

Next, we have 

Theorem 4. The nz(L y , d) are determined as follows. One has 

n z (L y ,d)=0 for d > L y , (4.2) 
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n z (L y ,L y ) = 1 



and 



n z (l,0) = l. 

All other numbers nz(L y , d) are then determined by the two recursion relations 

n z {L y + 1, 0) = n z (L y , 0) + n z (L y , 1) 

and 

n z (L y + l,d) = n z (L y ,d- 1) + 2n z (L y ,d) +n z (L y , d+ 1) 

for 1 < d < L y + 1 . 



(4.3) 
(4.4) 

(4.5) 



(4.6) 



Proof. The proof is similar to that for Theorem 2: we substitute for from eq. (2.1) in 
eq. (4.1). We obtain another equation by differentiating this with respect to q once; 
another by differentiating twice, and so forth up to L^-fold differentiations. This yields 
L y + 1 linear equations in the L y + 1 unknowns, n z (L y , d), d — 0, 1, L y . We solve this set 
of equations to get the n z (L y , d). □ 
A corollary is that 

n z (Ly,L y -l) = 2Ly-l . (4.7) 



From Theorems 3 and 4 we find a general formula for the numbers n z (L y , d): 

(T >x (2d+l) ( 2L„ 
n z {Ly,d) - 



(L y + d+ 1) VA, - d 



(4.8) 



for < d < L y , with n z (L y) d) = for d > L r For fixed d in the range < d < L y where 
n z (L y ,d) is nonvanishing, it has the leading asymptotic behavior 



n z (L y , d) ~ (2d + l^-^L-^' 2 A L y \l + 0(L- y l 



CIS -Z^ 7 



oc for fixed d 



(4.9) 



(The formal notation l + C^L^ 1 ) in (4.9) and in other asymptotic expansions is not intended 
to indicate the sign of the coefficient of the 0(L~ X ) term; here it is actually negative.) As a 
measure of the asymptotic behavior of n z (L y , d) when d, rather than being fixed, is a finite 
fraction of L y , we take the central value d = L y /2 and calculate that 



n z (L y , ^) ~ ^-1/2^1/242^+13-1(^+1) 1 + Q{L -u 

Note that for the special case d — 0, this reduces to 

n z (L y ,0)=C Ly 



oc 



(4.10) 



(4.11) 
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Table 3: Table of numbers nz{L y ,d) and their sums, Nz,g,\- See text for general formulas. 



L y I d -> 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Nz,Ly,\ 


1 


1 


1 




















2 


2 


2 


3 


1 


















6 


3 


5 


9 


5 


1 
















20 


4 


14 


28 


20 


7 


1 














70 


5 


42 


90 


75 


35 


9 


1 












252 


6 


132 


297 


275 


154 


54 


11 


1 










924 


7 


429 


1001 


1001 


637 


273 


77 


13 


1 








3432 


8 


1430 


3432 


3640 


2548 


1260 


440 


104 


15 


1 






12870 


9 


4862 


11934 


13260 


9996 


5508 


2244 


663 


135 


17 


1 




48620 


10 


16796 


41990 


48450 


38760 


23256 


10659 


3705 


950 


170 


19 


1 


184756 



where the Catalan number C n was defined in (3.19). In Table 3 we list the first few numbers 
nz{L y ,d) and the total sums N Zj l v ,x that will be calculated below in (5.2.2). 

Similar to Table 2, we can also make up a Catalan triangle, and shall find that Table 3 
is just the combination of two triangles. 

Certain (g-independent) relations between the riz(L y , d) can be derived by evaluating the 
sum (3.9) and its result (4.1) for q = and q = 1. Setting q = in this sums and using 
(2.17), we have 

Ly 

"£(-l) d n z (L y ,d)=0. (4.12) 

d=0 

Next we evaluate (3.9) and (4.1) for q — 1. Let L y = 3k + r with k > and r = 0, 1, or 2. 
Using (2.18), we have 

k k -i+e{r-D 
5>z(L„,3j)- n z (L y ,3j + 2) = l . (4.13) 

3=0 3=0 



We have obtained the following relations involving both nz(L y , d) and np(L y , d), as well 
as the total numbers Nz l a and Np l a 

n P (L y ,d) = J2(-iy ( Ly 7 l \n z {L y -j, d) (4.14) 

.7=0 
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n z (L y ,d) = J2 ■ ^npiLy-j^d) . (4.15) 

Since the total numbers of terms Np^ \ and Nz,l v ,x are sums of the np(L y , d) and nz(L y , d) 
given, respectively, by (3.2) and (3.3), it follows that 

NpM = E ( LV j ^ N 'M* ( 416 ) 



j"=o 

and 



^ 'L„ - 1 



= E \- ( 4 - 17 ) 

j'=o \ J / 

Note that for d = 0, our eq. (4.14) reduces to the relation (3.18) expressing the Motzkin 
number as a certain weighted sum of Catalan numbers, while eq. (4.15) reduces to the 
relation expressing the Catalan number as a weighted sum of Motzkin numbers (e.g., [67]) 

c »=e("; 1 )^- 1 =e("; 1 )^ < 4 - i8 > 

for n > 1 (with C = 1). Further, for the even and odd widths L y = 21 and L y — 2£ + 1, 

2 ±} (U - 1\ 

N z ,t,\ = E( -1 ) j ( ? W*-** ( 4 - 19 ) 
j=0 \ J J 

and 

2JV ZA a = E(- 1 )' ? N W-i+ ■ ( 4 - 2 °) 
i=o ^ J ' 



5 Determination of Np^ y ,\ and Nz,l v ,x for Cyclic and 
Mobius Strips of the Square and Triangular Lattices 
and Relation with Directed Lattice Animals 

In this section we shall use our results above to calculate the total number, Np t c,x, of Ap^j's 
in the chromatic polynomial (1.10) and the total number, Nz,g,x, of Az,g,/s in the full Potts 
model partition function for cyclic strips G of the square and triangular lattice. Since the 
individual numbers np(L y , d) were shown to be the same for the cyclic strips of the square and 
triangular lattices, and similarly for nz(L y , d), clearly it is also true that the respective total 



20 



numbers Np^,x and Nz,g,x are the same for cyclic strips of the square and triangular lattices. 
Furthermore, as a consequence of eqs. (3.4)- (3. 7), our use of cyclic strips to calculate the 
total numbers N p> g,x and N z ,g,x also yields these respective numbers for the Mobius strips 
of the square and triangular lattices. This is useful since, as will be seen in the next sections, 
the individual coefficients involve a larger set for Mobius strips of the square lattice, namely 
not just but also — c^; moreover, the coefficients for Mobius strips of the triangular 
lattice are not, in general, polynomials in q [42]. 



5-1 N P ^x 

For the total number of terms Np r \, from our theorems above we obtain the recursion 
relation 

Np, Ly+ i,x = 3N PtLytX - 2n P (L y ,0) . (5.1.1) 

We find that 

[— ] 

N P , Ly , x = 2S Ly = 2(L y -l)l± {j ^ L ~i ] m (5-1.2) 

where S n was given in (3.24). A generating function for these numbers is 

1 -I- 1/2 00 
Gn w {x) = ( T ^) " - 1 = £ N P , Ly , x x L « . (5.1.3) 

Ly = l 

From this, it follows that the number Np^ Ly ^\ grows exponentially fast with the width L y of 
the cyclic or Mobius strip of the square or triangular lattice, with the leading asymptotic 
behavior 

N P , Ly ,x ~ Ly 1 ' 2 3 L » as L y - oo . (5.1.4) 



5.2 N Z , G ,X 

A corollary of Theorem 4 is that 

N z , Ly+ i,x = 4iVz,L„A - 2n z (L y , 0) . (5.2.1) 

Summing the individual numbers nz(L y , d) to evaluate the total, (3.3), we have, for cyclic 
and Mobius strips of width L y of the square and triangular lattices, 



<k(L, + d+l)\L t -d) 



d=0 

2L y 

Ly 

As L y — > oo, Nz,l v ,x has the leading asymptotic behavior 



(5.2.2) 



z, Ly ,x ~ n-^L- 1 / 2 A L y [l + 0(L y x )] as L„ - oo . (5.2.3) 
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We give two comparisons of the asymptotic behavior of individual numbers nz{L y ,d) with 
the total, N z ,L y ,\- First, we compare the growth of riz{L y ,d) for fixed d (in the range 
< d < L y where nz(L y , d) is nonzero) with the growth of N z ,L y ,\'- 



n z (L y ,d) (2d+l) 



Z,Ly,X 



as L y — > oo for fixed d . 



(5.2.4) 



Second, we compare the relative growths of the central number nz(L y , -£) and the total, 



n z (L 



yi 2 



N Z,L y ,X 



/ 4 \L y + lr 

iw) [ l + 0{L ~y 

(0.769800.. .) Ly+1 
as L„ — > oo . 



1 + 



(5.2.5) 



5.3 Connection with Directed Lattice Animals 

For an arbitrary G, the Potts model partition function Z(G,q,v) gives information about 
certain graph-theoretic quantities describing G, as is clear from the Kasteleyn-Fortuin rep- 
resentation (1.5) and the equivalence in eqs. (12.1.6) and (12.1.9) to the Tutte polynomial 
and Whitney rank polynomial. An example of this information is illustrated by the zero- 
temperature antiferromagnetic special case P(G,q) = Z(G, q, — 1) in eq. (1.6), i.e., the 
chromatic polynomial, counting the number of proper vertex colorings of G. Other exam- 
ples are that the Tutte polynomials T(G,x,y) for various special values of its arguments, 
counts the number of spanning trees, spanning forests, connected spanning subgraphs, and 
spanning subgraphs of G (see appendix). It is also well known that the q — > 1 limit of the 
Potts model is related to bond percolation and the q — > 1 limit of a certain multisite Potts 
model is related to site percolation, as reviewed in [2]. 

Here we would like to report two new and very interesting relations with graph-theoretic 
quantities that we have found from our calculations. These involve directed lattice animals. 
To explain this, we recall that an animal A on a lattice or more generally a graph G is defined 
as a (finite) set of vertices in G with the property that any two vertices in A are connected 
by means of a path in G having all its vertices in A. Thus, any vertex in A is adjacent to 
another vertex in A. If the graph is a lattice, the animal is termed a lattice animal, and 
on an infinite lattice, the lattice animals are defined up to an overall lattice translation. A 
directed lattice animal A dir is an animal on a lattice A together with a special origin or 
root point O such that any vertex of A can be reached starting from O by an oriented path 
on A having all of its vertices in A. For example, a directed lattice animal on the square 
lattice could be defined by restricting the orientations of each step of the oriented path to be 
"eastward" or "northward", so that, starting at the origin, this animal would extend in the 
northeast quadrant and would not contain any backtracking steps going in the south or west 
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directions. (The animals that we consider here are also called site animals to distinguish them 
from bond animals; we shall leave this implicit.) Animals and directed animals are related, 
respectively, to percolation and directed percolation. Let us denote Na,a,ti an d N D a,a,u as 
the total number of lattice animals and directed lattice animals with n vertices, respectively, 
on the lattice A. For n — > oo, the numbers of lattice animals and directed lattice animals 
grow asymptotically like n~ e a n , where a depends on the lattice and, e.g., for 2D lattices, 
9 = 1/2 (e.g. [70]- [76] and references therein). Having given this background, we now state 
the relations that we have obtained. 

First, we find that the total number of distinct eigenvalues of the coloring matrix, i.e. the 
total number, Np t L \, of different Ap^'s in the chromatic polynomial (1.10) for cyclic strips 
with width L y of the square and triangular lattice (equal by (3.5) to the same number for 
the corresponding Mobius strips of the square and triangular lattices) is twice the number 
of directed lattice animals with n = L y vertices on the square (sq) lattice: 



We established this by observing that the generating function for directed lattice animals on 
the square lattice, which is known exactly [70, 71], is precisely (1/2) times the generating 
function (5.1.3). 

Second, we find that the analogous total number, Nz t L y ,x, of different eigenvalues \z,g,/s 
appearing in the full Potts model partition function (1.9) for these cyclic and Mobius strips 
with width L y of the square and triangular lattices is twice the number of directed lattice 
animals with n = L y vertices on the triangular lattice: 



and using our calculation (5.2.2). It should be emphasized that, as we have shown above, the 
numbers N PtLyjX and N z ,l v ,\ on the respective left-hand sides of (5.3.1) and (5.3.2) apply for 
cyclic and Mobius strips of both the square and triangular lattice, whereas the right-hand 
side of (5.3.1) is specific to the square lattice and the right-hand side of (5.3.2) is specific to 
the triangular lattice. 

Parenthetically, we mention some other relations. The enumeration of directed lattice 
animals on the square and triangular lattices was shown in [71] to be connected with the 
hard-square lattice gas model [72]. As discussed in [75, 76], no exact result has been obtained, 
analogous to those in [70, 71], for the number of directed lattice animals on the honeycomb 
lattice. Directed lattice animals are connected with directed percolation, and the latter has 
been related to a kind of chiral Potts model [77]. Compact lattice animals are the subset 
of lattice animals with no unoccupied interior sites; directed compact lattice animals have 
been related to the q — > oo limit of the Potts model and have been shown to have numbers 
that grow less rapidly than the leading a n growth in the numbers of regular directed lattice 
animals [78]. 



Np t L y ,\ — 2N D A,sq,L y ■ 



(5.3.1) 



Nz,L y ,\ — 2N D A,tri,L y ■ 

We established this by recalling the known result [70, 71] 



(5.3.2) 




(5.3.3) 
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6 Determination of np^Mb(L y , d, ±) for Mobius Strips of 
the Square Lattice 

Mobius strips differ in their global topology from cyclic strips. For the Mobius strips of 
the square lattice that have been studied, the coefficients are polynomials in q, but the 
polynomials that occur for the strip of width L y arise from the set ±c^, where < d < L y /2 
for even L y and < d < (L y + l)/2 for odd L y , rather than the set d d \ < d < L y as 
in the case of the cyclic strips. In passing, we note that for Mobius strips of the triangular 
lattice, we have found, originally for width L y = 2 in [42] and later for width L y = 3 in [46], 
that in the case of chromatic polynomials, some of the coefficients are not polynomials, but 
rather, algebraic functions of q [42, 46]; similarly, some of the coefficients in the full Potts 
model partition function for the Mobius strip of the triangular lattice are not polynomials 
in q [47]. This is related to the fact, discussed before [46, 47], that the Mobius strip of 
the triangular lattice involves a seam, i.e., it is not translationally homogeneous in the 
longitudinal direction, whereas the Mobius strip of the square lattice does not have a seam 
and is translationally homogeneous in the longitudinal direction. Hence, for Mobius strips, 
we only consider the square lattice here. 

Since both signs of the occur, it is necessary to define a larger set of numbers for 
Mobius strips of the square lattice. Let np t Mb{L y ,d,±) denote the number of terms X/s in 
the expression (1.10) for P(sq, L y x L x , FBC y , TPBC X , q) with coefficients cp t L y ,Mb,j = 
respectively, where Mb denotes Mobius. The notation cp^ y ,Mb,j reflects the fact that these 
coefficients depend on L y and the Mobius boundary conditions, but are independent of L x . 
The sums and differences for each d are defined as 

np,Mb,tot(L y , d) = n Pt Mb(L y , d, +) + n PjM b(L y , d, -) (6.1) 

and 

An PtMb (L y , d) = n Pt Mb(L y , d, +) - n PtMb (L y , d, -) . (6.2) 
We first obtain a general theorem for the sum of the coefficients: 

Theorem 5. 



"^' A |0 for even L y 

Cp, Ly ,Mb = }^ °P:Ly,Mb,j = 2^ An PiM b(L y ,d)c [ = < p(T Ly+1 q ) for odd L y ( 6 - 3 ) 



3=1 d=0 

where P(T n , q) was given in (1.13) and 



_ . for even L y 

dmax " < for odd L y ■ [bA) 



Proof. Our method for proving this theorem is inspired by coloring matrix methods [43]. 
The Mobius strip involves a reversed-orientation periodic longitudinal boundary condition. 
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We can think of constructing such a strip by cutting a cyclic strip, reversing the orientation 
of one of the ends, and gluing these ends together again. For a strip with L y even, let us 
label the vertices on the two ends as 1, 2, .., L y \ then the Mobius boundary condition means 
identifying vertex 1 with vertex L y , vertex 2 with vertex L y — 1, and so forth. As regards 
the coloring constraint, in the case L y = 2, this effectively produces a subgraph consisting 
of a single vertex with an edge forming a loop. In the case L y = 4, this produces a vertex 
1 = 4 with two edges connecting to a vertex 2 = 3, which in turn has a loop connected to 
it, and so forth for higher even values of L y . Because of the loop that appears in each 
even-L y case, the chromatic polynomial for coloring this subgraph vanishes identically. 
(The value of d max given in (6.4) follows from this graphical construction.) This proves the 
theorem for the case of even L y . For odd L y , consider first the case L y = 3; here the 
Mobius boundary condition identifies vertex 1 with vertex 3 and leaves vertex 2 invariant. 
Hence, as regards the coloring, it effectively produces a subgraph consisting of the vertex 
1 = 3 connected by two edges with the vertex 2. The chromatic polynomial for the coloring 
of this subgraph is not sensitive to the multiple edges and hence is P(T2, q) = q(q — 1). 
Again, this subgraph construction yields the value of d max for odd L y in (6.4). In a similar 
manner, for higher odd values of L y , the Mobius boundary condition leads to a subgraph 
with (L y + l)/2 vertices forming a chain, with each interior pair connected to the next by 
two edges. The chromatic polynomial for the coloring of this subgraph is 
P(T( Ly+1 y 2 , q) = q{q — l)^ -1 )/ 2 . This completes the proof of the theorem. □ 

Two corollaries of this theorem are as follows. 
Corollary 1 

n P ,Mb(L y , d, +) = n Pt Mb(L y , d, -) for even L y (6.5) 

Proof. If L y is even, then in order for the terms of the highest power in q to cancel so as to 
yield a sum of as in (6.3), it is necessary that n Pj Mb(L y , d max , +) = n Pj Mb(L y , d max , — ), 
where d max was given above in (6.4). But then in turn, in order for the terms of degree 
d max — 1 to sum to zero, it is necessary that 

n P> Mb{L y , d max - 1, +) = n P)M b{L y , d max - 1, -), and so forth for all powers. □ 
Corollary 2 

An PtMb (L y ,d) = n P { Ly ^ - ,d) for odd L y > 3 (6.6) 

Proof. This follows from eq. (6.3) by the same kind of argument that was used to obtain 
the n P (L y ,d). Given that the sum C Pt L y ,Mb satisfies (6.3), this uniquely determines the 
An Pj Mb(L y , d)'s just as (3.10) determined the n P (L y ,d). □ 
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In order to determine the numbers of coefficients c^ d \ i.e., np : Mb(Ly, d, ±), for a Mobius 
strip of width L y , we start with the chromatic polynomial for the cyclic L y strip and 
determine how the coefficients change when one changes the longitudinal boundary 
condition from cyclic to Mobius. The next two theorems determines this (the proofs are 
given after the second of these two theorems). 

Theorem 6 Consider an L y x L x strip of the square lattice with (FBC yi PBC X ). As before, 
the coefficients cp t c,j are made up from the set with < d < L y . When one changes the 
longitudinal boundary condition from cyclic to Mobius, the following respective changes of 
coefficients of even degree d = 2k and of odd degree 2k + 1 occur: 



.(0) 



,(2fc) 



±C 



(fc-1) 



±c<°> 



,1 < k < 



2 \ 



c (2*+i) _> ±c (*+i) ,0<A;< 



(6.7) 
(6.8) 
(6.9) 



where in eqs. (6.8) and (6.9), we again use the notation [u] to denote the integral part of v. 
Thus, if there are n P (L y , 0) coefficients cp jGj - of the terms (\j) m of the form in eq. 
(1.10) for P(sq(L y x L x ; FBC y , PBC x ),q), then the respective coefficients multiplying the 
terms {\j) m in the chromatic polynomial for the corresponding Mobius strip, 
P(sq(L y x L x ; FBC y , TPBC X ), q), are either of the form or — c^°\ and so forth for the 

with d > 0, as specified by eqs. (6.8) and (6.9). 

Theorem 7. With the same premise as in Theorem 6, we have 



1 1 

np t Mb{L y , 0, ±) = -np t Mb,tot(Ly,0) = - 



np(L y ,0) + np(L y ,2) for L y even (6.10) 



np,Mb(L y , d, ±) 



2 n P,Mb,tot(L y ,d) — - 



n P (L y , 2d - 1) + n P (L y , 2d + 2) 



for L y even and 1 < d < 



rip t Mb(L y ,d,±) = for L y even and d > 



n PtMb (L y , 0, ±) = i n P (L y , 0) + n P (L y , 2) ± n P ( Ly ~^ 1 , 0) 



Ly + 2 



for L y odd 



(6.11) 

(6.12) 
(6.13) 



np,Mb(L y , d, ±) 



n P (L y , 2d-l) + n P (L y , 2d + 2)± n P ( Ly + - , d) 
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for L y odd and 1 < d < — ^ - 



np t Mb(Ly, d, ±) = for L y odd and d > 



L y + 3 



(6.14) 
(6.15) 



Proofs. We prove these theorems using an inductive argument. From the chromatic 
polynomials for the cyclic and Mobius strips of the square lattice with width L y = 2 [18], 
one knows the transformation rule for the coefficients with < d < 2 as one goes from 
cyclic to Mobius longitudinal boundary conditions, namely, c*- -* — > ±c^°\ — > ±c^\ and 
c*- 2 - 1 — > ic*- - 1 , where the sign information will not be needed. For L y = 3, we know the 
numbers n P (L y , d) for the cyclic strip from the previous theorem. We also know the 
differences of the numbers of coefficients of each degree for the L y = 3 Mobius strip, 
Anp t Mb(3,d); these are determined by the relation (6.3). In particular, we have 
Anp M6 (3, 0) = n P (2, 0) = 1, Anp M fc(3, 1) = n P (2, 1) = 2, and Anp M b(3, 2) = n P (2, 2) = 1. 
From the transformation rules obtained so far, it follows that the sources for coefficients of 
degree for the L y = 3 Mobius strip are the coefficients of degree and 2 for the 
corresponding cyclic strip. Adding these, we have n Pt Mb,tot(3, 0) = np(3, 0) + n P (3, 2) = 5. 
This relation for the sum of n PjMb (3, 0, +) and npMft(3, 0, — ), together with the relation 
AnpMb(3, 0) = 1 for their difference, enables us to solve for each of these numbers, and we 
get n P: Mb(3, 0, +) = 3, n P: Mb(3, 0, — ) = 2. From the transformation rule c 1 ^ — > ic^ found 
previously, we infer that for the Mobius strip, n PtMbjtot (3, 1) = n P (3, 1) = 4. Combining this 
with the relation AnpMb(3, 1) = 2, we solve to get n Pt Mb(3, 1, +) = 3 and 
npMfe(3, 1, — ) = 1. We next consider npMft(3, 2, ±). We know that AnpMb(3, 2) = 1, and, 
given the transformation rules obtained so far, there is only one source for terms with 
degree coefficients of degree 2 in the chromatic polynomial for the L y = 2 Mobius strip, 
namely the single term with a coefficient. Hence, n Pt Mb,tot(3,2) = 1, and so, knowing 
the sum and difference, we get n P< Mb(3, 2, +) = 1 and n Pt Mb(3, 2, — ) = 0. Since there are no 
sources for any higher-degree coefficients, we have npMj>(3, d, ±) = for d > 3. Thus, we 
have both determined all of the numbers n P) Mb{3, d, ±) and the transformation rule for the 
coefficients of next higher degree, namely — > ±c^. 

Proceeding to L y = 4, we start with the knowledge of the np(4, d) for the cyclic strip and 
of the differences AnpMb(4, d) = 0. For coefficients of degree in the chromatic polynomial 
for the Mobius strip, the transformation rules indicate two sources, namely the coefficients 
with degree or 2 in the corresponding cyclic strip. We shall show, it a posteriori, that 
these are the only sources, so we infer that n P: Mb,tot(^, 0) = np(4, 0) + np(4, 2) = 12. 
Combining this with the relation AnpMb(4, 0) = 0, we obtain 

^p,Mfe(4,0, ±) = (l/2)np 5 Mfe,tot(4, 0) = 6. Next, for the coefficients of degree 1, we start with 
those from the cyclic strip, np(4, 1) = 9. But we know that there must be another source 
because this is an odd number and, by itself, would give the unacceptable non-integral 
result for npMb(4,0, ±) = (l/2)n P; Mb,tot{4=, 0). The only possible source is the coefficients of 
degree 4 in the cyclic strip, so we infer the next transformation rule: c( 4 ) — > ±c < - 1 \ This 
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yields np t Mb,tot(^, 1) = np(4, 1) + np(A, 4) = 10 and hence np t Mb(^, 0, ±) = 5. Finally, for the 
coefficients of degree 2 in the chromatic polynomial for the L y = 4 Mobius strip, the source 
is the coefficients of degree 3 in the cyclic L y = 4 strip, so np t Mb,tat(^, 2) = np(4, 3) = 4, 
whence np 5 Mb(4, 2, ±) = 2. Since there are no sources for any higher-degree coefficients, we 
have np i Mb(4, d, ±) = for d > 3. Thus again we have determined all of the numbers 
n P,Mfe(4, <i, ±) and also the next higher transformation rule for c*- 4 -*. 

It should now be clear how one proceeds iteratively: at each higher width L y , one uses the 
previously established transformation rules, and obtains the transformation rule for the 
next higher degree coefficient, to determine all of the numbers np t Mb{L y , d, ±). This 
completes the proof of Theorems 5 and 6. □. 

Two corollaries of these theorems are 

n P:Mb (L y , y,±) = y foreven L y ( 6 - 16 ) 

n PM b(L y , Ly ^ ~ , +) = 1 , n Pt Mb(L y , Ly ^ 1 , -) = for odd L y (6.17) 

Values of the numbers np t Mb(L y , d, ±) for the first several widths, 2 < L y < 10, are given in 
Table 4. 

Further corollaries involve (g-independent) relations between the n P (L y ,d). As before 
for the cyclic strips, we derive these by evaluating (6.3) for q — 0,1, and 2. For example, 
setting q = and using (2.17), we have (for L y > 2 where Mobius strips are defined) 

dmax f Lt, r T 

^ I ror even 

2^(-l) d An P: Mb(L y ,d) = where d max =< $ +1) y (6.18) 

d =o { ~^2~ tor ocld L v 

It is straightforward to derive similar relations among the Anp t Mb(L y , d) by evaluating (6.3) 
for q — 1 and q = 2. 



Another result pertains to the detailed structure of the chromatic polynomials for the 
cyclic, as compared with Mobius strips of a given width L y , as the length L x gets large. As 
discussed in [44], because the cyclic and Mobius strips of a given width have the same number 
of vertices and edges, the coefficients of the leading powers of q are the same. We recall this 
result. For m greater than a minimal value 4 , the cyclic and Mobius strips of a given (G s ) m 
have the same number of vertices n, edges (bonds) e, girth g (length of minimum closed 
circuit on G s ) and number k g of circuits of length g. One has n = t s m where t s depends on 
G s . Writing 

n— 1 

P((G a ),q) = Y,(-l) i hn- i <r i (6-19) 

3=0 

4 For example, for m < 4 (m < 3) the L y = 2 cyclic (Mobius) square strips exhibit special behavior 
regarding g and k g : for m larger than these respective values, they both have g = 4 and k g = m. 
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Table 4: Table of numbers np,Mb(L y , d, ±) for Mobius strips of width L y . For each L y value, the entries in 
the first and second lines are np,Mb(L y , d, +) and np y Mb(L y , d,—), respectively. Blank entries are zero. The 
last entry for each value of L y is the total Np^ y ,Mb,\- 



L y i (d,+) - 


o,+ 


1,+ 


2,+ 


3,+ 


4,+ 


5,+ 




(d,-) - 


o,- 


1,- 


2,- 


3,- 


4,- 


5,- 




2 


1 
1 


1 
1 










4 


3 


3 
2 


3 
1 


1 








10 


4 


6 
6 


5 
5 


2 
2 








26 


5 


16 
14 


15 
11 


8 
5 


1 






70 


6 


38 


35 


20 


3 










38 


35 


20 


3 






192 


7 


100 


100 


64 


15 


1 








96 


91 


56 


11 






534 


8 


256 


263 


178 


49 


4 








256 


263 


178 


49 


4 




1500 


9 


681 


739 


535 


178 


24 


1 






672 


718 


514 


165 


19 




4246 


10 


1805 


2027 


1539 


574 


96 


5 






1805 


2027 


1539 


574 


96 


5 


12092 
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and using the results that [13, 20] h n -j = for < j < g — 1 (whence h n — 1 and /i n _i = e) 

and = (^Ij — fc 9 , it follows that for m greater than the above-mentioned minimal 

value, these h/s are the same for the cyclic and Mobius strips of each type G s . For a given 
G s , asm increases, the h n ^/s for the cyclic and Mobius strips become equal for larger j. 
Our relations (6.7)-(6.9) make clear which coefficients c Gs ( Ly )j in (1.10), and hence which 
coefficients hc s (L y ),j in (6.19), remain equal for the cyclic versus Mobius strips. 

7 Determination of the numbers nz,Mb(L yi d, ±) for Mobius 
Strips of the Square Lattice 

In this section we generalize the calculations of the previous section to the full Potts model 
partition function on Mobius strips of the square lattice. (The restriction to the square lattice 
will be implicit in the notation.) Let nz,Mb(L y , d, ±) denote the number of terms \z,L y ,j in 
the expression (1.9) for Z(L y x L x , FBC y ,TPBC x , q,v) with coefficients cz,L y ,Mb,j = 
respectively, where Mb denotes Mobius. The notation cz,L y ,Mb,j reflects the fact that these 
coefficients are independent of L x . The sums and differences for each d are defined as 



nz,Mb,tot(L y , d) — Tiz,Mb(L y , d, +) + nz,Mb(L y , d, — ) 



and 



Anz,Mb(Ly, d) — nz,Mb(L y , d, +) — n z ,Mb(L y , d, — ) 
We first have the two theorems 

Theorem 8. 



(7.1) 
(7.2) 



Nz,Ly,\ 



d m ax 



Cz,L y ,Mb = ^2 c z,L y ,Mb,j — ^2 An z ,Mb(Ly, d)c 

j=l d=0 



(d) 



A/2 



for even L„ 



q (L y +i)/2 for odd Ly 



(7.3) 



Theorem 9. For the Mobius strips of the square lattice, 

n z ,Mb(L v ,0,±) = \ nz(L y ,0) + n z (L y ,2) ±n z (^,0)] for even (7.4) 



n z,Mb{Ly, d, ±) 



n z (L y , 2d-l)+ n z (L y , 2d + 2)± nz(y , d) 
for L„ even and 1 < d < — 



n z ,Mb( L y, 0, ±) = ^ n z (L y , 0) + n z (L y , 2) ± n z ( Ly ^ 1 , 0) 



for L y odd 



(7.5) 
(7.6) 
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1 



L„ + l 



n ZMb (L y ,d,±) = - n z (L y , 2d- 1) + n z (L y ,2d + 2) ±n z ( — - — , d) 



for L y odd and 1 < d < . (7.7) 

nzMLy, d,±) = for d > [^^1 (7.8) 



The proofs of these theorems are analogous to those for the analogous Theorems 5 and 7 
for chromatic polynomials for Mobius strips of the square lattice. 

Two corollaries of these theorems are 

n z ,Mb(L y , = L y , nz,Mb( L y, -) = L y ~ 1 for even L y (7.9) 

n z ,Mb(L y , Ly ^ - , +) = 1 , n ZtMb (L y , Ly ^ - , -) = for odd L y (7.10) 

Values of the numbers n Zt Mb(L y , d, ±) for the first several widths, 2 < L y < 10, are given in 
Table 5. 

Two relations among the n Zt Mb(L y , d) can be derived as before, by evaluating the sum 
(7.3) for q = and g = 1. For example, from the evaluation for q = we have (for L y > 2) 

g(-l) d (An z , M6 (L„d)) = where ^ = ^ °^^ w (7.11) 
It is straightforward to obtain a corresponding relation via the evaluation of (7.3) for q = 1. 

8 Some Examples 

Since some of the notation above is complicated, it is worthwhile to illustrate our general 
calculations with specific exact solutions. These will also be useful in a discussion of the 
special cases of the infinite-temperature limit of the Potts model and the zero-temperature 
limit of the Potts ferromagnet below. 

8.1 Chromatic Polynomials 

For the L y = 3 cyclic strip of the square lattice, iVpL =3 / \ = 10 and, with the abbreviation 
sq3FP for sq, L y = 3, FBC y , PBC X , the chromatic polynomial is [42] 

10 

P(sq,3 x L x ,FBC y ,PBC x ,q) = J2 c ^FpA^ q3 FP,j) Lx (8.1.1) 

3=1 
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Table 5: Table of numbers nz,Mb(L y , d, ±) for Mobius strips of width L y . For each L y value, the entries in 
the first and second lines are nz,Mb(L y , d, +) and nz,Mb{L v , d, respectively. Blank entries are zero. The 
last entry for each value of L y is the total Nz,L v ,Mb,\- 



L y i (d,+) - 


o,+ 


1,+ 


2,+ 


3,+ 


4,+ 


5,+ 




(d,-) - 


o,- 


1,- 


2,- 


3,- 


4,- 


5,- 


Nz,L y ,Mb,\ 


2 


2 


2 














1 


1 










6 


3 


6 


6 


1 












4 


3 










20 


4 


18 


16 


4 












16 


13 


3 








70 


5 


61 


54 


20 


1 










56 


45 


15 








252 


6 


206 


180 


80 


6 










201 


171 


75 


5 






924 


7 


722 


651 


335 


42 


1 








708 


623 


315 


35 






3432 


8 


2542 


2360 


1336 


224 


8 








2528 


2332 


1316 


217 


7 




12870 


9 


9082 


8766 


5367 


1148 


72 


1 






9040 


8676 


5292 


1113 


63 




48620 


10 


32644 


32668 


21270 


5432 


480 


10 






32602 


32578 


21195 


5397 


471 


9 


184756 
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where 



A 



sq3FP,l 



-1 



A S(? FP,(6,7) - 



1/2' 



(8.1.2) 

(8.1.3) 
(8.1.4) 
(8.1.5) 
(8.1.6) 

(8.1.7) 



Kq3FP,2 — q — 1 
^sq3FP,3 — Q ~ 2 

A sg 3FP,4 = q — 4 

A sg 3FP,5 = — (q — 2) 2 

(g - 2)(q 2 - 3g + 5) ± {(g 2 - 5g + 7)(g 4 - 5g 3 + llg 2 - 12g + 8)} 
and X sq FP,j, j = 8, 9, 10, are the roots of the cubic equation 

£ 3 + (2g 2 - 9q + 12)^ 2 + (g 4 - 10g 3 + 36g 2 - 56g + 31)^ 

-(g - l)(g 4 - 9g 3 + 29g 2 - 40g + 22) = . 

The coefficients are 

C S g3PP,l = C (3) 

c sq 3FP,j = c (2) for j = 2,3,4 
c S g3PP,j = c (1) for j = 5, 8, 9, 10 
c sq 3FP,j = c (0) = 1 for j = 6, 7 . 

For the numbers n P (L y , d), one sees that n P (3, 0) = 2, np(3, 1) = 4, n P (3, 2) = 3, n P (3, 3) = 
1, and the total number Np t L y ,\ = 10, in agreement with Table 1. The exact solution for the 
chromatic polynomial for the L y = 3 cyclic strip of the triangular lattice, given in [46], has 
the same N P)Ly ^ and np(L y , <f)'s, although eight of the ten Ap^/s are different. 

As an illustration of our results for chromatic polynomials for Mobius strips of the square 
lattice, we again take the L y = 3 case. The exact solution for the chromatic polynomial (with 
the abbreviation sg3M6 for sq,L y = 3, FBC y ,TPBC x ) is [44] 



(8.1.8) 

(8.1.9) 
(8.1.10) 
(8.1.11) 
(8.1.12) 
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P(sq,3 x L x , FBC y ,TPBC x ,q) = ^ cp jSq3M b,j(Xp,sg3Mb,j) Lx 

3=1 

where \p, sq 3Mb,j = Xp,sq3FP,j V j, in accordance with the general result (3.4) and 

(2) 

Cp, S q3Mb,l — C 

cp, S q3Mb,j = c (1) for j = 8, 9, 10 



(8.1.13) 



Cp <S q3Mb,5 — 



c P,sq3Mb,j 



= c (0) for j = 3, 6, 7 



(8.1.14) 

(8.1.15) 
(8.1.16) 
(8.1.17) 
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cp,s q 3Mb,j = -c (0) for j = 2,4 . (8.1.18) 

Thus, for the numbers np(L y , d, ±), we have np(3, 0, +) = 3, np(3, 0, — ) = 2, n P (3, 1, +) = 3, 
np(3, 1,— ) = 1, np(3,2,+) = 1, and n P (3,2,— ) = 0, with np(3, d, ±) = for d > 3, in 
agreement with Table 4. 

The fact that the coefficients in the chromatic polynomial for the Mobius strip of the 
triangular lattice are not, in general, polynomials in q, was shown via the exact solution for 
the lowest width, L y = 2 in [42] (with L x = m): 



P(tri, 2 x L x , FBCy, TPBC X , q) = -1 + [(q - 2) 2 } m 

[(m-l)/2] 



(g-l)(g-3) 
V9^4^ 



(Aj2Mb,3) m — (A t 2M6,4) r 



-l + (g-2) 2m -2 1 - m (g-l)(g-3) g x j (5 - 2,)— (9 - 4^ (8.1.19) 



where 



A t2 Mfe,(3,4) = 2 5 - 2g ± a/9 - 4q J . 



(8.1.20) 



As we have discussed in [46], this can be attributed to the feature that the Mobius strip of the 
triangular lattice has a seam, and hence is not translationally homogeneous in the longitu- 
dinal direction, in contrast to the Mobius strip of the square lattice, which is translationally 
homogeneous in the longitudinal direction. 

8.2 Potts Model Partition Function 

As an illustration for the full Potts model partition function, we recall the calculation of this 
function for the L y = 2 cyclic strip of the square lattice. Here, for Nz,l v ,x, we have [15, 16] 
N z ,2,\ = 6, and (with the shorthand sq2FP for sq, L y = 2, FBC y , PBC X ) 



where 



A 



Z,sq2FP,(3A) - 



Z(sq, 2 x L x , FBC y , PBC X , q,v) = J2 Cz,s q 2FP,j(\z ,sq2FP,j) 

3=1 



\ -2 

AZ,sq2FP,l — V 

>>Z,sq2FP,2 = V{V + q) 

q + v(v + 4) ± (v 4 + 4v 3 + 12v 2 - 2qv 2 + 4qv + q 2 ) 1/2 



^Z,sq2FP,(5fi) — -j{Ts\2 ± >/ Rsi2 ) 



with 



T S i2 = v 3 + 4v 2 + 3qv + q 2 



R S12 = v ° + 4v° - 2qv" - 2q 2 v 3 + \2v" + 16^ d + 13g z w z + 6q% + g 4 . 



2„.2 



(8.2.1) 

(8.2.2) 
(8.2.3) 

(8.2.4) 
(8.2.5) 

(8.2.6) 
(8.2.7) 
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and 

cz, SQ 2FP,i = c {2) (8.2.8) 

cz,s q 2FP,j = c (1) for j = 2, 3, 4 (8.2.9) 

cz^fpj = c (0) for j = 5,6 . (8.2.10) 

Thus, for the individual numbers, n^(2, 0) = 2, n^(2, 1) = 3, n z (2, 2) = 1, and n z (2, d) = 
for d > 3, in agreement with Table 3. 

As an example of our general results for the structure of the Potts model partition function 
on Mobius strips of the square lattice, we display the exact solution for width L y = 2 [15, 16] 
(with the abbreviation sq2Mb for sq,L y = 2, FBC y ,TPBC x ): 

6 

Z(sq,2 x L x ,FBC y ,TPBC x ,q,v) = ^ c ZySq2Mbj (X ZySq2Mb ^ (8.2.11) 

where \z, sq 2Mb,j = ^z,sq2FP,j V j, in accord with the result (3.6) and 

c z , sq 2Mb,j = c (1) for j = 3,4 (8.2.12) 

cz, sg2 MM = -c {1) (8.2.13) 
cz^Mbj = c (0) for j = 5,6 (8.2.14) 
c M m M = -c {0) (8.2.15) 

Thus, n Z)Mb (2,0, +) = 2, n ZjM fe(2,0, -) = 1, n Z;M6 (2, 1,+) = 2, and n ZjMb (2, 1,+) = 1, in 
agreement with Table 5. 



8.3 Other Limits 

In the text above we have given results for structural properties of the general Potts model 
partition function Z(G, q, v) and for the special case of the zero-temperature Potts antiferro- 
magnet, Z(G,q,v = —1) = P(G,q). Another special case is infinite-temperature, i.e., v = 0. 
In this case, for any graph G, the Potts model partition function reduces to Z(G, g, 0) = q n ( G \ 
as indicated above in (1.8). For the strip graphs under consideration here, the equality in 
(1.8) corresponds to the property that at v — all but one of the A^gj's in eq. (1.9) vanish, 
and the nonvanishing \z,G,j necessarily has the value q Ly and the coefficient = 1, so that 
the right-hand side of (1.9) becomes q L y Lx = q n . One can see this explicitly in the exact 
solution given above for the Potts model partition function for the L y = 2 strip of the square 
lattice. 

Secondly, one may consider the zero-temperature limit of the Potts ferro magnet, v — > oo. 
This is quite different from the T = limit of the Potts antiferromagnet; among other 
differences, the ground state entropy per site is zero for the ferromagnet rather than being 
nonzero as it is for the antiferromagnet with sufficiently large q. In this limit T — > 0, the 
partition function of the Potts ferromagnet grows like a e ^ where a was defined in (1.3) and 
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e(G) was also defined in the introduction as the number of edges (bonds) in G. Hence, it is 
convenient to use the reduced partition function Z r , defined by 

Z r (G, q, v) = a~ e(G) Z{G, q, v) = u e(G) Z{G, q, v) (8.3.1) 

which has the finite limit Z r — > 1 as T — > 0. For the cyclic or Mobius strip graphs of the 
lattice A of interest here we thus write 

N Z ,G,\ 

Z r (A,FBC y ,(T)PBC x ,q,v) = u<^ £ c z ,gAXz,gj) l ' 

3=1 

= c z,gAXz,gj,u) L ° (8.3.2) 

where as before we use the shorthand notation G — A, FBC y , (T)PBC X , and 

\z,G,j,u = « e(0)/£ "A z , Gj - • (8.3.3) 

For example, for cyclic or Mobius strips of the square lattice, e = (2L y — 1)L X , so the 
prefactor in (8.3.3) is u 2Ly ~ l , while for the corresponding strips of the triangular lattice, 
e = (3L y — 2)L X , so the prefactor is u 3Ly ~ 2 . For the L y — 1 circuit graph, 

Az,c,i,u = 1 - w , A ZiC)2 ,n = 1 + (g - l)u . (8.3.4) 

so that both of these two Az,gj,u' s remain important in the T = limit. In contrast, for the 
L y = 2 cyclic or Mobius strips of the square and triangular lattices, our exact solutions in 
[16, 47] show that, of the six Xz,G,j,u s , two are finite and nonzero for u — > while the other 
four vanish. For example, for the cyclic/Mobius strip of the square lattice, in the vicinity of 
the zero-temperature point u — 0, one has [16] 

Xz,sq2FP,l,u = U-2U 2 + U 3 (8.3.5) 

\z, sq2 FP,2,u = u + (q- 2)u 2 + (1 - q)u 3 (8.3.6) 
and the Taylor series expansions 

\z, sq 2FP,3,u =l~u 2 + 2(q- 2)u 3 + 0( M 4 ) (8.3.7) 

\z, sq 2FP A ,u = u+(q- A)u 2 + (7 - 3q)u 3 + 0(u 4 ) (8.3.8) 

\z, sg 2FP, 5 ,u = 1 + (q - l)u 2 [l + Au + 0(m 2 )] (8.3.9) 

\z, sq 2FP,e,u = u + 2(q- 2)u 2 + (q 2 -7q + 7)u 3 + 0(u 4 ) . (8.3.10) 

Hence, 

, \z, sq 2FP, 3 ,u _ , X z , aq 2FP,j _ [ 1 for j = 5 f , 

AZ, sq 2FP,3,u AZ,s q 2FP,3 L U 101 — Z ' ^' 

In this sense, one can define an effective, reduced, Nz,g,x, and this has the value 2, somewhat 
analogous to the value Np t o,\ = 4 for the T = limit of the Potts antiferromagnet (1.10); 
however, a difference is that when one reinserts the prefactor in (8.3.2) to get back the actual 
partition function, all of the six \z,g/ s do contribute. 
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9 Coefficients for Other Strip Graphs 



9.1 Strips with Torus or Klein Bottle Boundary Conditions 

Like the cyclic and Mobius strips, the strips with torus or Klein bottle boundary conditions 
have periodic or twisted periodic longitudinal boundary conditions, respectively. We first 
give two theorems for the sum of the coefficients for these types of strip graphs (which are 
the analogues of the theorems yielding eqs. (3.10) and (6.3) for the cyclic and Mobius strips). 
We have 

Theorem 10. The sum of the coefficients in (1.16) for a strip of the square or triangular 
lattice of width L y and (PBC y , PBC X ) = torus boundary conditions is (independent of L x ) 

Np } L y ,torus,\ 

Cp j L y ,torus= ^ Cp t L y ,torus,j — P(CL y ,q) (9.1.1) 
3=1 

where the chromatic polynomial of the circuit graph C n is 

P(C n ,q) = (q- l) n + (q - 1)(-1)» . (9.1.2) 

Proof. This follows by noting that the transverse slice of the strip graphs of the square 
and triangular lattices of width L y is the circuit graph Cl v , and the number of ways of 
proper colorings of this graph is given by the chromatic polynomial in (9.1.2). □ 

Theorem 11. The sum of the coefficients in (1.16) for a strip of the square or triangular 
lattice of width L y is (independent of L x ) and (PBC y ,TPBC x ) = Klein bottle (KB) 
boundary conditions is 

Np,L y ,KB,\ 

Cp, Ly ,KB= Yl C P,L y ,KB,j = (9.1.3) 
3=1 

Proof. This follows by observing that this sum corresponds to the proper coloring of an 
effective graph (defined as in our study of Mobius strips) formed from the transverse slice 
of the strip, and in this case the effective graph always involves at least one loop, i.e. an 
edge that connects a vertex to itself (specifically, the effective graph involves one loop if L y 
is odd and two loops if L y is even). But the chromatic polynomial of a graph with a loop 
vanishes identically. □ 

Theorems 10 and 11 generalize the observations made for these sums of coefficients for 
the exact solutions of the chromatic polynomials of the L y = 3 torus and Klein bottle graphs 
of the square lattice in [45] and [46] to arbitrary L y values where these graphs are defined, 
i.e. L y > 3. 
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On the basis of the exact solutions for the chromatic polynomials for the width L y = 
3 strips of the square [45] and triangular [46] lattices with (PBC y , PBC X ) = torus and 
(PBC y ,TPBC x ) = Klein bottle boundary conditions, we can make several further remarks. 
First, since it was found that N PjG: x = 8 for G = sq(3 x L X) PBC y) PBC X ) but the different 
value, Np y c,\ = 11 for tri(3xL x , PBC y , PBC X ), these strips behave fundamentally differently 
than the cyclic and Mobius strips of these two lattices, for which the value of Np t a,\ was 
the same for a given L y . Secondly, the exact solutions in [45, 46] show that the coefficients 
that enter in the chromatic polynomial for the L y = 3 strips of the square and triangular 
lattice with torus or Klein bottle boundary conditions are polynomials in q but are not of 
the form = L^^)- Indeed, in contrast to the situation for cyclic and Mobius square- 
lattice strips and cyclic triangular-lattice strips, the coefficient of degree d is not unique, 
i.e., there can be more than one type of coefficient of a given degree d in q. For example, 
for the L y = 3 square-lattice strip with torus boundary conditions, the coefficients that are 
quadratic functions of q are q(q — 3), (l/2)q(q — 3), (q — l)(q — 2), and (l/2)(g — l)(q — 2). 
Another difference, again illustrated by the chromatic polynomial for the L y = 3 strips of 
the square and triangular lattice is that the coefficients can have zeros outside the range 
< q < 4, in contrast to the c^ d \ For example, one term in the chromatic polynomial for 
the L y = 3 strip with toroidal boundary conditions is (g 3 — 6q 2 + 8g — 1)(— l) Lx for the 
square lattice [45] and (l/3)(g - l)(q 2 - 5q + 3)(-2) Lx for the triangular lattice [46]; the first 
coefficient has a zero at q ~ 4.11 and the other coefficient has a zero at q ~ 4.30. Given 
that the exact solutions for the chromatic polynomials on torus graphs in [45, 46] show that 
the coefficients, which are multiplicities of eigenvalues of the coloring matrix, are not of the 
form for strip graphs of the square and triangular lattices, it is not clear to us how to 
relate this with the remark in [29] that certain related eigenvalue multiplicities d r for torus 
boundary conditions are of the form sin(r#/2)/ sin(6/p) where p is 1 (2) if r is even (odd) 
and 9 satisfies (2.25). 

The exact solutions for the chromatic polynomial in [45] for the square lattice and in [46] 
for the triangular lattice showed that if one starts with torus boundary conditions, cuts the 
tube graph and reglues the ends with reversed orientation to form Klein bottle boundary 
conditions, the number of Ap^./s does change (is reduced): 



Thus, at least for this width, if one uses Klein bottle boundary conditions, then the value of 
Np,G,\ is the same, namely five, for both square and triangular strips. As is evident from the 
exact solutions [45, 46], the five coefficients are the same, although the Xp t G,/ s are different. 
These coefficients are 




(9.1.4) 



(9.1.5) 



(9.1.6) 



{1, q-1, -(q-1), 



(g-l)(g-2) q(q-3) 
2 ' 2 



} 



(9.1.7) 
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Since the coloring matrix applies directly for the lattice strip graphs with torus boundary 
conditions, eq. (1.18) applies for the chromatic polynomials of these graphs. 

9.2 Strips with Free Longitudinal Boundary Conditions 

Previously, exact solutions for chromatic polynomials on strip graphs of regular lattices 
with (FBCy, FBC X ) = open boundary conditions were given in [33, 40]; in these cases, as 
discussed before, the coefficients cpc-j are not, in general, polynomials in q. This can be 
seen immediately from eqs. (2.14) or (2.19) of [40]. An analogous statement holds for the 
full Potts model partition functions on these open strips [15, 16, 47]. This property can be 
understood as a consequence of the fact that the chromatic polynomial is not the trace of 
an m'th power of the coloring matrix for these strips, but rather is given by 

P(G,q) = (c\T m \c') (9.2.1) 

where c and d denote colorings of the ends of the strip. A similar statement holds for the 
full Potts partition function Z(G,q,v) in terms of the corresponding matrix 7z. 

As discussed before [36, 41, 42, 44, 15, 46], for a strip graph of some lattice with given 
transverse boundary conditions, the term \p,G,j that is dominant in the physical region of q 
(including sufficiently positive integers and denoted as region Ri in our previous work [14]) 
is independent of the longitudinal boundary conditions. (The same applies for the Xz,gj m 
Z(G,q,v) that is dominant in the physical paramagnetic phase [15, 16, 47].) In particular, 
this means that the A p t G,j that is dominant in region Ri is the same for the cyclic (or Mobius) 
and open strip graphs of a given lattice type A, thereby establishing a certain connection 
between these cyclic (or Mobius) and open strip graphs. Similarly, it means that the Xp,G,j 
that is dominant in region R 1 is the same for the torus and Klein bottle strips of this lattice 
A on the one hand, and the cylindrical strip graphs on the other. Now it was also shown [15] 
that for graphs with periodic or twisted periodic boundary conditions, the coefficient of this 
dominant term is unity. If this term is a root r^ s of an (irreducible) algebraic equation of 
degree dt, then the theorem on symmetric functions of roots of an algebraic equation implies 
that for the cyclic strips, these roots rg, s appear in P(G, q) in the form of the sum YltLi r eh 
i.e., they all have the coefficient unity [44]. This sets an upper bound on dg(A,L y ), namely 

np(L y ,0)>d e (A,L y ) (9.2.2) 

where we have included the dependent of di on the lattice type A and the strip width, and 
we recall from eq. (3.4) that dg(A, L y ) is the same for the cyclic and Mobius strips of a given 
lattice A. For cyclic and Mobius strips of the square lattice, dg(sq, 1) = 1, dg(sq,2) = 1 
[18], de(sq : 3) = 2 [42], and de(sq,A) = 3 [48], while for the cyclic and Mobius strips of the 
triangular lattice, de(tri,2) = 1 [42] and de(tri,3) = 2, de(tri,A) = 4 [46]. This information 
is summarized in tables given in [46, 48] along with related results. Our calculations above 
yield, for the relevant values of np(L y , 0) (see Table 1): np(l, 0) = np(2, 0) = 1, np(3, 0) = 2, 
and rap (4, 0) = 4. Hence, evidently, the inequality (9.2.2) is realized as an equality for the 
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triangular strips and the square strips for L y = 1,2,3, but as a strict inequality for the 
L y = 4 square strips. Further, for the cases that we studied, we found that 

de(A, L y ) C y C , )M b = de(A, L y ) open (9.2.3) 

where dg(A, L y ) open = de(A, FBC y , FBC X , L y ). When we obtained the results in the present 
paper for np(L y ,d), we noticed the interesting connection that for the cyclic and Mobius 
strips of the triangular lattice, for which (9.2.2) holds as an equality, it is also true that one 
has the fourfold equality 

n P (L y , 0) = de(tri, L y ) cyc ^ M b. = dg(tri, L y ) open = N P ^ ri ( Ly ^ open ^ (9.2.4) 

i.e., the number of Ap^'s with coefficient c<°) = 1 in the chromatic polynomial for the cyclic 
or Mobius strips of the triangular lattice, which was equal to the degree of the algebraic 
equation one of whose roots was the dominant Xp,G,j in was a l so equal to the total 
number of terms in the chromatic polynomial for the open strip of the triangular lattice 
with the same width. We confirmed eq. (9.2.4) for widths L y = 2,3,4, based on the exact 
solutions in [33, 42, 46]. In terms of the generating function formalism of [33, 34, 40], the 
denominator V of the generating function, as a polynomial in the auxiliary variable, has 
the maximal degree, np(L y ,0), allowed by (9.2.2) for these open strips of the triangular 
lattice but does not, in general, for the open strips of the square lattice. As we noted in 
[46, 48] , this difference in behavior is related to the fact that the Mobius strips of the square 
and triangular lattices are different; in particular, the Mobius strip of the triangular lattice 
involves algebraic nonpolynomial coefficients. Further, even for the strips of the square 
lattice, where the first two terms of (9.2.4) are not, in general, equal (as our solutions for 
chromatic polynomials for L y = 4 showed), we still found 

de(sq, L y ) cyc , t Mb. = N PjSq ( Ly -) j0penj \ ■ (9.2.5) 

That is, we find that the equality (9.2.5) holds for cyclic or Mobius and open strips of the 
square lattice for 1 < L y < 4, as well as for the corresponding strips of the triangular lattice 
with 2 < L y < 4; in particular, it holds even when np(L y , 0) is not equal to (is greater than) 
di(A,L y ). For reference, for open strips of the triangular lattice with L y > 2, one can infer 
from calculations in [29] that N Zjtri(Ly)j0penjX = C Ly ; furthermore, Np jtri(Lyhopent x = M Ly _ x 
[49] . Another generalization would be that for triangular strips, the value of the degree in 
eq. (9.2.3) satisfies di(tri, L y ) cyc ^ Mb = di(tri, L y ) open = Ml i. This would require that, in 
terms of the generating function [33], the denominator V does not factorize. 

Previously, coloring matrix methods have been applied to the honeycomb (he) lattice to 
derive rigorous bounds on W [23, 24]. In the present context of strips, we observe that, 
in addition to the elementary result Np t hc,L y =2,o P en,x = de(hc,L y = 2) open = 1, our exact 
solutions have yielded the results 

Np t hc,L y =3,open\ = de(hc, L y = 3) open = 3 (9.2.6) 

from [33] and 

Afp t hc,L y =2,cyc.,Mb.,x = de(hc, L y = 2) cyC ;Mb = 1 (9.2.7) 
from [41]. In view of eq. (9.2.6), for this strip of the honeycomb lattice, rip(L y = 3, 0) ^ 
d e (hc, L y = 3) open , in contrast to the case with the square and triangular lattice strips. 
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9.3 Strips with Cylindrical Boundary Conditions 

Exact solutions for chromatic polynomials for cylindrical strips, i.e. with boundary con- 
ditions of the form (PBC y , FBC X ) of the square and triangular lattices were presented in 
[34, 36] , and, again, the coefficients are, in general, not polynomials in q, for the same reason 
as in the case of the open strips. And again, the same statement holds for the full Potts 
model partition function. (Recent calculations of chromatic polynomials for wider strips 
with open and cylindrical boundary conditions include [46, 48, 49].) 



10 Some General Geometric Identities 

In this section we present some useful identities between Potts model partition functions on 
different types of lattice strips that follow from basic geometrical considerations. We have 
applied these as checks in our previous calculations. Consider a family of strip graphs of the 
square or triangular lattice with fixed width L y and arbitrary length L x with some set of 
longitudinal (x) and transverse (y) boundary conditions, denoted as BC = (BC y , BC X ). By 
the transfer matrix argument given in [16] , it follows that the Potts model partition function 
for this strip is of the form 

Nz,A,L y ,BC 

Z(A,L y x L x ,BC,q,v) = ]T c z , AtLutB cj(Xz,A,Ly,BCj) L ' ■ (10.1) 

i=i 

We recall [16] that the total number of terms, A r z,A,L y ,BC, the coefficients cz,A,Ly,BC,j, and 
the \z,K,L y ,BC,i s are independent of L x . In our previous works, we have studied the L x — > oo 
behavior of these Potts model partition functions. However, clearly, for fixed L y and L x , by 
switching the x and y axes, one can equivalently view the strip as the length-L^ member of 
a family of strip graphs with fixed width L x and length L y . b This yields a number of useful 
identities. We proceed to describe these 



10.1 Strips with (FBC y , FBC X ) 

For strips of the square and triangular lattices with (FBC y , FBC X ) (open) boundary condi- 
tions, the interchange of the x and y axes leaves the boundary conditions the same, and one 
gets the following identity (where we write F y F x = (FBC y , FBC X ) to save space) 

N Z ,A,L y ,F y ,F x N Z ,A,L x ,F x ,F y 

j(^Z,A,L x ,F x ,F y ,j) V ■ 

3=1 3=1 

(10.1.1) 



5 Note that for other lattices such simple relations do not, in general, hold. For example, if one considers 
strips of the honeycomb lattice, constructed as a brick lattice, then rotating a open horizontal L y x L x strip 
of bricks by 90°, one does not get an equivalent L x x L y strip of horizontally oriented bricks. 
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Taking T = for the Potts antiferromagnet, one obtains the corresponding identity for 
chromatic polynomials 



Np,h,Ly,F y ,F x Np tA ,L x ,Fy,F x 

^ Cp,A,Ly,F y ,F x ,j(^P,A,Ly,Fy,F x ,j) Lx = ^ Cp^ Lx ,F x ,F y ,j(^P,A,L x ,F x ,F y ,j) Ly . 
3=1 j=l 

(10.1.2) 

Using the results in [33] , one can construct various illustrations of this. One of the simplest is 
to take the open strip of the square lattice with width L y = 2, for which Np tSq ^^FBC y ,FBC x — 1 



and P(sq,2 x L x , FBC y , FBC X , q) = q(q - l)(q 2 - 3q + 3) 



L x -1 



For L x = 3, this chro- 



matic polynomial must be equal to that for the width L y = 3 open strip (which has 
Np tSq fi,FBCy,FBC x = 2) with length L x = 2, for which (using [33] and eq. (2.15) of [40]) 
for general L x , 



P(sq,3x L x , FBCy, FBC y , q) 



(^sq3FF,l — A S g3i?F,2. 
\L x — 2 



■X 



(A sq 3 FF fl\ sq3FF: i + A S g3i?i7 ) i)(A sg 3i? j p i i)^ a: 2 — (A sq3 pp :0 \ sq 3Fp t 2 + A sq3 pF : i)(\ sq3 FF j2 ) Lx 2 

(10.1.3) 



where 



and 



A sq3F F,o = q(q - 1)(<? 2 - 3g + 3) 2 
A sq3FF ,i = -q(q - l) 3 (g 3 - 6g 2 + 13q - 11) 



(10.1.4) 
(10.1.5) 



A 



sq3FF,(l,2) 



(q-2)(q 2 -3q + 5)± (q 2 - 5q + 7)(g 4 - 5<f + llq 2 - 12q + 8) . (10.1.6) 
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10.2 Cyclic and Cylindrical Strips 

We have denoted strips of fixed width L y and arbitrary length L x vertices with (FBC y , PBC X ) 
and (PBC yi FBC X ) boundary conditions as cyclic and cylindrical, respectively. Topologi- 
cally, these are the same; the distinction was made because we were, in particular, interested 
in the large L x — > oo limit, which, for cyclic strips means that the length of a circuit goes 
to infinity, while for the cylindrical strips, the length of the circumference of the cylinder is 
fixed while its length goes to infinity. Clearly, a cyclic strip of the square or triangular lattice 
with a fixed width L y and length L x is the same as a cylindrical strip of the given lattice 
with width L x and length L y . Letting A denote the type of lattice, square or triangular as 
before, and writing 

Nz,A,Ly,Fy,P X 

Z(A, L y x /„.. FBC y , PBC X ) = £ c ZA l^,p^(A w ^,p^) L * (10.2.1) 
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we have the identity 

Nz,A,L y ,F y ,P x N Z ,A,L x ,P y ,F x 

^ CZ,A,Ly,F y ,P x ,j(^Z,A,Ly,F y ,P x ,j) Lx = ^2 c Z,A,L x ,P y ,F x ,j(>^Z,A,L x ,P y ,F x ,j) Ly ■ 
3=1 3=1 

(10.2.2) 

Again, for the special case of the T = Potts antiferro magnet, we have the resultant identity 
for the chromatic polynomial 

N Pi A,Ly,Fy,P X Np iA ,L X ,Py,F X 

^ c P,A,Ly,Fy,P x ,j(^P,A,Ly,F y ,P x ,j) Lx = ^ Cp , A , L x , Py , F x , j ( A p A , L x , Py , F x , j ) L V ■ 
3=1 3=1 

(10.2.3) 

Perhaps the simplest illustration of this general identity is the relation between the chromatic 
polynomial for the width L y = 2 cyclic strip of the square lattice and the width L y = 3 
cylindrical strip of this lattice. For the L y = 2 cyclic strip, one has Np tSq ^xL x ,F y ,p x = 4 and 

P(sq,2xL x ,FBC y ,PBC x ,q) = (q 2 -3q+l) + (q-l) [(3-g) L * + (1 -g) L *] +(g 2 -3g + 3) L * . 

(10.2.4) 

For the L y = 3 cylindrical strip, one has 

P(sq, 3 x L x , PBC y , FBC X , q) = q(q - l)(q - 2)(q 3 - 6q 2 + Uq - 13) L *- L . (10.2.5) 
Thus, the identity (10.2.3) yields 

(g 2 -3g+l) + (g-l)[(3-g) 3 + (l-g) 3 ] +(g 2 -3g + 3) 3 = 

q(q - l)(q - 2)(q 3 - 6q 2 + Uq - 13) . (10.2.6) 

While this illustration is quite simple, more complicated cases involve identities between 
powers of algebraic roots of different types. 

10.3 Strips with Torus Boundary Conditions 

For a strip graph of the square or triangular lattice with (PBC y , PBC X ) = (P y , P x ) (torus) 
boundary conditions, the identities read 

Nz,A,Ly,Py,P x Nz,A,L x ,P Xl P y 

^ CZ,A,Ly,Py,P x ,j(^Z,A,Ly,P y ,P x ,j) Lx = ^ C Z , A,L X ,P X ,P y ,j ( ^Z, A,L X ,P X ,P y , j ) Ly 
3=1 3=1 

(10.3.1) 

N P: A,Ly,Py,P X ^ Pl A,L Xl P X ,Py 

^ CP,A,Ly,Py,P X ,j(^P,A,Ly,Py,P X ,j) LX = ^ ^P , A,L X ,P X ,Py ,j ( ^P, A,L X ,P X ,Py ,j ) ^ 

3=1 3=1 



10.3.2) 
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The simplest illustration of this is for the family of strips of the square lattice with width L y = 
2, for which the torus degenerates to a cyclic strip, and one has P{sq, 2xL x , PBC y , PBC X , q) = 
P(sq, 2 x L x , FBCy, PBC X , q), given in (10.2.4). Let us next consider the strip of the square 
lattice with width L y = 3 and length L x = 2. This graph is related by the identity to two 
other strip graphs: (i) the cyclic strip of the square lattice with L y = 2 and L x = 3 (see 
eq. (10.2.4)) and (ii) the equivalent cylindrical strip of the square lattice with L y = 3 and 
L x = 2 (see eq. (10.2.5)). Indeed, when one evaluates the exact solution given in [45] for 
L x = 2, one finds that it is equal to the right-hand side of eq. (10.2.6). 

11 Conclusions 

In this paper, using a general formula, in terms of Chebyshev polynomials of the second 
kind, for the coefficients that occur in the g-state Potts model partition function on cyclic 
strips of the square and triangular lattices and on Mobius strips of the square lattice, to- 
gether with general formulas for sums of these coefficients, we have determined the num- 
ber, nz(L y ,d), of \z,g,/s with coefficient in Z(G,q,v) and the total number, Nz,g,x, 
for these strips. Results are also given for the analogous numbers np(L y ,d) and Np t L y ,\ 
for the zero-temperature Potts antiferromagnet partition functions, i.e., chromatic polyno- 
mials for these strips. The results for the total numbers, Nz,l v ,x and Np t L y ,\, apply for 
both cyclic and Mobius strips of both the square and triangular lattices. Among other 
connections, we found that np(L y ,0) = np(L y — 1,1) = Ml v -i, the Motzkin number; 
nz(L y ,0) = Cl v , the Catalan number; the exact expression for Np t L y ,x in eq. (5.1.2); and 
the relations N P)Ly ^ = 2N DAySqyLy ; and N z ,l v ,\ = 2N DA:trijLy , where N DA)A>n denotes the 
number of directed lattice animals on the lattice A. We also found the asymptotic growths 
Nz,L v ,\ ~ L y 1 A Ly and N P)Lyi \ ~ L y ' 3 y as L y — > oo. Some commens about other 
lattice strips were made. In addition, we presented some useful general geometric identities 
for Potts model partition functions. 
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12 Appendix 

12.1 Connection Between Potts Model Partition Function and 
Tutte Polynomial 

The Potts model partition function Z(G,q,v) is related to the Tutte polynomial T(G,x,y) 
as follows. The graph G has vertex set V and edge set E, denoted G — (V, E). A spanning 
subgraph G' is defined as a subgraph that has the same vertex set and a subset of the edge 
set: G' = (V,E') with E' C E. The Tutte polynomial of G, T(G,x,y), is then given by 
[7]- [9] 

T(G,x,y)= ^2{x-l) k ^- k{G \y-lf G ^ (12.1.1) 

G'CG 

where k(G'), e(G"), and n{G') = n(G) denote the number of components, edges, and vertices 
of G', and 

c(G') = e(G') + k(G') - n(G') (12.1.2) 

is the number of independent circuits in G 1 (sometimes called the co-rank of G'). Note that 
the first factor can also be written as (x — l) r ( G ) _r ( G ') ; where 

r(G) = n(G) - k{G) (12.1.3) 

is called the rank of G. The graphs G that we consider here are connected, so that k(G) = 1. 
Now let 

x = l + - (12.1.4) 

v 

and 

y = a = v + l (12.1.5) 
so that q — (x — l)(y — 1) = (x — l)v. Then 

Z(G,q,v) = (x-l) k ^ G \y-ir^T(G,x,y) . (12.1.6) 
There is also a connection with the Whitney rank polynomial, R(G,£,r)), defined as [4, 10] 

R(G,t, V )= £f (G V (G,) (12.1.7) 

G'CG 

where the sum is again over spanning subgraphs G" of G. Then 

T(G, x, y) = {x- l) r ^R(G, £ = (x - rj = y - 1) (12.1.8) 

and 

Z(G,q,v) = q n ^R(G,£ = ^ V = v) . (12.1.9) 
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Note that the chromatic polynomial is a special case of the Tutte polynomial: 

P{G, q) = q HG)^k(G)+n(G) T ( G ^ x = 1 _ q ^ y = ^ (12.1.10) 

(recall eq. (1.6)). 

For a recursive family of graphs, such as the strip graphs considered in this paper, com- 
prised of m repetitions of a basic subgraph, the Tutte polynomial has the form [16] 

N T,G,X 

T(G,x,y)= '■i,,,i\ l ,;,r (12.1.11) 

where, from (12.1.6), one has the relation 

\r,Gj=v- L y\ z ,Gj (12-1.12) 

so that 

N t ,g,x = N Z)G)X . (12.1.13) 
It is convenient to extract a common factor from the coefficients: 

= f^f • (12.1.14) 

Of course, although the individual terms contributing to the Tutte polynomial are thus 
rational functions of x rather than polynomials in x, the full Tutte polynomial is a polynomial 
in both x and y. Given the relation (12.1.6), if one defines fiT{L y , d) as the number of terms 
\r,L y ,j in T(L y ,q,v) that have as their reduced coefficient cr,L y ,j = c^ d \ then these are the 
same numbers: 

n T (L y ,d) = n z (L y ,d) . (12.1.15) 

Thus, Tables 3 and 5 apply equally to the structure of the Tutte polynomials for the cyclic 
and Mobius strips of the square and triangular lattices. 

For a given graph G = (V,E), at certain special values of the arguments x and y, the 
Tutte polynomial T(G,x,y) yields quantities of basic graph-theoretic interest [9]-[ll]. We 
recall some definitions: a spanning subgraph was defined at the beginning of the paper; a 
tree is a connected graph with no cycles; a forest is a graph containing one or more trees; 
and a spanning tree is a spanning subgraph that is a tree. We recall that the graphs G that 
we consider are connected. Then the number of spanning trees of G, N S t(G), is 

N ST (G)=T(G,1,1) , (12.1.16) 
the number of spanning forests of G, N S f(G), is 

N SF {G)=T(G,2,1) , (12.1.17) 
the number of connected spanning subgraphs of G, Ncssg(G), is 

N CSSG (G) = T(G,1,2) , (12.1.18) 
and the number of spanning subgraphs of G, Nssg(G), is 

N SSG (G)=T(G,2,2) . (12.1.19) 
These connections have been used in [79], [80]. 
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12.2 Determinants of Coloring Matrices 

In this section we list the determinants of coloring matrices for some lattice strips. For the 
chromatic polynomial this is 

det T P (G) = JjAp, Gj - = J] &P,aj) e "' aj ■ ( 12 - 2A ) 

3=1 3=1 

For the Potts model partition function the analogous determinant is 

JVz n z,g,x 

detT z (G) = J]A z>0j - = J] ( A *Ai) c * lOJ • ( 12 - 2 - 2 ) 
j=i j"=i 

Let us define the shorthand notation 

Dp[G) = det Tp(G) , D Z (G) = detT z (G) . (12.2.3) 

For the products of eigenvalues contributing to (1.10) for the chromatic polynomials and 
(1.9) for the Potts model partition function for cyclic and Mobius strips of the square lattice 
we find 

D P (sq, L y = 1, FBCy, PBC X ) = D P ({C}, q) = (-l) 9 " 1 ^ - 1) (12.2.4) 

(where {C} refers to the circuit graph) For the L y = 2 cyclic and Mobius strips of the square 
lattice, from [18], 

D P (sq,L y = 2,FBC y ,PBC x )= [(3 - q)(l - g)] "~V - 3q + 3) (12.2.5) 

and 

D P (sq,L y = 2,FBC y ,TPBC x ) = (3 - ^(l - ^"^V - 3g + 3) . (12.2.6) 
For the L y = 3 cyclic and Mobius strips of the square lattice, from [42], 

Dp(sq, L y = 3,FBC y ,PBC x ) = (_i)(«-2)( 9 2 -3 9+ D ((? _ i)(«-i) 2 ( g _ 2 )« a -*- 1 x 

(g - 4) g2 - 3,?+1 (g 3 - 6g 2 + 13g - ll)(g 4 - 9g 3 + 29g 2 - 40g + 22) 9 " 1 (12.2.7) 

and 

D P (sq, L y = 3,FBC y ,TPBC x ) = {-lf-^ +2 {q - \y~\q - 2)-( 2 *" 3 > x 

(9 - 4) _1 (g 3 - 6g 2 + 13g - ll)(g 4 - 9g 3 + 29g 2 - 40g + 22) 9 " 1 . (12.2.8) 

The corresponding determinants for L y = 4 can be obtained from [48] but are quite lengthy, 
so we do not list them here. 
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For the Dz determinants we have 

D z (sq, L y = 1, FBCy, PBC X ) = D Z ({C}, q) = v q ~\q + v) (12.2.9) 

and, from [16], 

D z (sq, L y = 2, FBCy, PBC X ) = v 2qiq ~ l \v + l) q (v + q) 2q (12.2.10) 

and 

D z (sq, L y = 2, FBC y ,TPBC x ) = v 2{q ~ l \v + l) q {v + qf . (12.2.11) 

Note that, except for the case of the circuit graph, sq, L y = 1, PBC X , for the value v — — 1, 
where the Potts model partition function reduces to the chromatic polynomial, Z[G,q,v = 
— 1) = P(G,q), some eigenvalues in the product contributing to Z vanish, and hence this 
product vanishes. If one extracts these vanishing eigenvalues, then, of course, the rest yield 
the same product as for P(G, q). 

For cyclic strips of the triangular lattice we have, from [42] 

D P (tri, L y = 2, FBCy, PBC X ) = (q - 2) 2q (12.2.12) 

and from [46, 47], 

D P (tri, L y = 3, FBCy, PBC X ) = (_i)(4-2)( 9 2 -3 9+ i) ((? _ 2 )<?( 2 <?- 1 )( g - 3)i(i- 1 ) (12.2.13) 

D z (tri,L y = 2,FBC y ,PBC x ) = v 2qiq - l \v + l) 2q {v + q) 2q . (12.2.14) 

The determinant D P {tri,L y = 4, FBC y , PBC X ) can be calculated from our exact solution 
in [46]; however, it is sufficiently lengthy that we do not include it here. The analogous 
determinants for the Mobius strips of the triangular lattice can also be calculated from the 
exact solutions that we have given [42, 46]; however, they are more complicated, since some 
coefficients are algebraic, rather than polynomial, functions of q. 

For the strip graphs of the square and triangular lattices with (PBC y , PBC X ) = torus 
and (PBCy, TPBC X ) = Klein bottle boundary conditions, we have, from the exact solutions 
in [45] and [46], 

D P (sq,L y = 3,PBC y ,PBC x ) = (-l) q3 - 6q2+llq -\q - l) i3 =^ 1 x 

( 5 _ 2) q2+q ~\q - A)( q '^ q - 2 \q - 5)^(g 2 - 7q + l^X 

(g 3 - 6q 2 + Uq - 13) (12.2.15) 

and 

D P (sq,L y = 3,PBCy,TPBC x ) = (q - 1)'^^^ (q - 5)^ (q 2 - 7q + l^x 

(q 3 - 6q 2 + Uq - 13) . (12.2.16) 
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12.3 Some Other Coloring Matrix Results for Lattice Strip Graphs 



Since the trace of a coloring matrix vanishes (recall (1.18)), the sum of eigenvalues, each 
multiplied by its multiplicity also vanishes, as indicated in eq. (1.19), for the cyclic and torus 
strips, for which (1.14) applies directly. Using coloring methods, we obtain the following 
general formulas for sums for other types of strips. These agree with our previous exact 
solutions in [42, 46, 48]. 



Np, aq {Ly),Mb,\ 



E = {f- 1){t - 3 « + 3 > T ~ (12.3.1) 

j=i < y 



N P,tri(L y ),Mb,\ 



Cp,tri(Ly),Mb,j^P,tri{L y ),Mb,j ~ V L y > 2 . (12.3.2) 

3=1 

The equations also apply to the corresponding strips with Klein bottle (KB) boundary 
conditions, with Mb — > KB, and L y > 3. The polynomial q 2 — 3q + 3 in (12.3.1) is D 4 in 
our previous notation, where 

D k ^(-l) s ( k ~ 1 )q k - 2 - s . (12.3.3) 
s=0 V s / 

From the exact solutions for the chromatic polynomials in [42], we have, for the cyclic strip 
of the kagome lattice consisting of a succession of hexagons with two interleaved triangles 
per hexagon of arbitrary length (denoted as having width L y = 2 in [42]) we have 

6 

CpmjXpmj = q(q ~ 1) 2 (<7 - 2) 2 . (12.3.4) 

j'=i 

Additional families of graphs are provided by the homeomorphic expansions of the cyclic 
and Mobius L y = 2 strips of the square lattice, in which one adds k — 2 degree-2 vertices 
to the upper and lower horizontal edges of each square, for k > 3. These may be viewed 
as strips of p-sided polygons, with p = 2k, with each successive pair of polygons sharing 
one edge. For k = 2 and k — 3, these families are the L y = 2 cyclic and Mobius strips of 
the square and honeycomb lattice, respectively. (In the latter case, the honeycomb lattice 
is represented as a brick lattice with the bricks oriented horizontally.) Following [41], we 
denote the cyclic (eye) and Mobius (Mb) strips of this type, of length L x = m as (CK) k m cyc 
and (Ch)k, m ,Mb, where Ch denotes "chain". Exact solutions for the chromatic polynomials 
of these families were given in [41]. From the construction of these graphs, it is clear that 
Ej c P,(Ch) k , cyc , 3 = q(q - 1) and J2j c P,(Ch) k , Mb}j = hold, as in eqs. (3.10) and (6.3). For the 
other sums (with X(ch),k,cyc,j = \ch),k,Mbj), w e have 

4 

E ^Chh^Ach)^ =q 2 -3q+l + 2{-l) k+ \q - l)D k+l + D 2k (12.3.5) 
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where Dk was defined in (12.3.3), and 

4 

E c (^,mmW,mm = -l + 2(-l) fe (g-l) J D fc + J D 2fc . (12.3.6) 
j'=i 

For k = 2, the sum (12.3.5) is zero, as in (1.18), and the sum (12.3.6) has the value q(q — 1) 
as in the L y = 2 special case of (12.3.1). As an example of the values for higher-p polygonal 
strips, for k — 3, i.e., the strip of the honeycomb lattice, (12.3.5) has the value q(q — l) 3 and 
(12.3.6) has the value q(q - l)(q - 2) 2 . 



For the full Potts model partition functions, from our exact solutions in [16, 47] we find 

6 

^ CZ,sq2,cyc,j^Z,s Q 2,cyc,j = q(v + if (v + q) (12.3.7) 
3=1 



^ Cz,sq2,Mb,j^z,sq2,Mb,j = q(v 3 + 3v 2 + 3v + q) (12.3.8) 

3=1 

6 

^ cz,tri2, C j/ C ,iAz,tri2, CJ/c ,i = g(u + 1) 2 (^ 2 + 2v + g) (12.3.9) 
i=i 

6 

^ C Z ,tri2,Mb,j^Z,tri2,Mb,j = q(v + 1)(^ 3 + 3^ 2 + 3v + q) . (12.3.10) 
3=1 

As is evident, for the special case v = — 1 where the Potts model partition function reduces 
to the chromatic polynomial, these equations reduce to their analogues for the respective 
chromatic polynomials. 

In terms of the Tutte polynomials, these formulas read 

6 

^2 CT,sq2,cyc,j^T,sq2,cyc,j = Xtj 2 (12.3.11) 
3=1 

6 

C T ,sq2,Mb,j^T,sq2,Mb,j = X + JJ + TJ 2 (12.3.12) 

3=1 

(which is the Tutte polynomial for the graph known as the "thick link" with three edges, 
the planar dual to the circuit graph C 3 ) 

6 

c T,tri2,cyc,j^T,tri2,cyc,j — ( x + U)y 2 (12.3.13) 

3=1 
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and 

6 

^2 C T,tri2,Mb,j>^T,tri2,Mb,j = y(x + X) + y 2 ) ■ (12.3.14) 
3=1 
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